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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3910 


THE RESPONSE OF AN AIRPLANE TO RANDOM 
ATMOSPHERIC DISTURBANCES 1 
By Franklin W. Diederich^ 

SUMMARY 


The statistical approach to the gust-load problem "which consists in 
considering flight through turbulent air to be a stationary random process 
is extended by inclu ding the effect of lateral variation of the instan- 
taneous gust intensity on the aerodynamic forces. The forces obtained in 
this manner are used in dynamic analyses of rigid and flexible airplanes 
free to move vertically, in pitch, and in roll. The effect of the inter- 
action of longitudinal, normal, and lateral gusts on the wing stresses 
is also considered. 

The method of analyzing the rigid-body motions is similar to that 
used for analyses of the dynamic stability of airplanes, in that the 
equations of motion are referred to stability axes and expressed in terms 
of conventional stability derivatives. The method of analyzing the 
dynamic effects of structural flexibility consists in an extension of a 
numerical-integration approach to the static aeroelastlc problem and is 
In a form which offers the possibility of calculating divergence and 
flutter speeds with relatively little additional effort. 

The mean-square values, correlation functions, and power spectra of 
seme of the aerodynamic forces required in this type of analysis are cal- 
culated for one special correlation function of the atmospheric turbu- 
lence. It Is shown, for instance, that if the span Is relatively large 
compared with the integral scale of turbulence, the mean-square lift and 
root bending moment directly due to the gust are substantially reduced 
when the differences in instantaneous intensity of the turbulence along 
the span are taken into account. However, if the motions of the airplane 
are taken Into account the mean-square root bending moment may be Increased 
as a result of these differences. Also, the mean-square pitching moment 
is shown to be substantially increased if the tail length is relatively 

^This report represents, except for some minor changes, a thesis 
submitted in May 195^ in partial fulfillment of the requirements for 
the degree of Doctor of Philosophy at the California Institute of 
Technology, Pasadena, Calif. 

%ow at the Langley Aeronautical Laboratory, Langley Field, Va. 
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large compared with the scale of turbulence . Finally, the wing stresses 
due to longitudinal, normal, and lateral gusts are shown to be statis- 
tically independent under certain conditions. 


INTRODUCTION 


The local velocity fluctuations acting on an airplane flying through 
atmospheric turbulence are functions of time defined only in a statistical 
sense and, hence, constitute a stochastic or random process. Consequently, 
the responses of the airplane, whether they are motions (linear or angular 
displacement, velocities, or accelerations), forces (lift, pitching moment, 
bending moment, and so on), stresses, or any other phenomena determined 
by the turbulence, can also be known as functions of time in only a sta- 
tistical sense. 

This paper is concerned with the statistical characteristics of those 
responses which have a bearing on the loads and stresses experienced by 
the airplane; although other problems such as those relating to passenger 
comfort or to the stability of the airplane as a gun platform can be 
treated in the same manner, they will not be considered here. 

The first approaches to the gust-load problem which use the statis- 
tical techniques developed for stationary random processes appear to be 
those of references 1, 2, and 3. (An earlier investigation concerned 
with the motions of an airplane in turbulent air is reported in ref. 4.) 

The fundamentals of these approaches are discussed in some detail in 
reference 1, and mention is made therein of Investigations in other fields 
of engineering and physics that have dealt with the problem of deducing 
the statistical characteristics of the output or response of a dynamic 
system from those of its input. The mean-square normal acceleration of 
a rigid airplane free to move in one degree of freedom, namely, vertical 
motion, is calculated in reference 2. In addition to being rigid, the 
airplane is implicitly assumed to be small enough for all its components 
to experience the same gust velocity at any instant of time. This means 
that the span of the airplane must be small compared with the integral 
scale of atmospheric turbulence, which on the basis of the available 
knowledge concerning the properties of the atmosphere (ref. b, for instance) 
appears to be in the order of several hundred to 1,000 or 2,000 feet; 
that is, the span of the airplane must be less than about 100 feet. 

The purpose of the present paper is. to extend this approach to large 
flexible airplanes free to move in all directions. As used herein, the 
terms "small" and "large" airplane refer 'to airplanes which are very small 
and not very small, respectively, compared with the integral scale of 
turbulence; thus, an airplane flying in a wide variety of atmospheric 
conditions may be "small" under certain conditions and "large" -under 
others. Similarly, the terms "rigid" airplane and "flexible" airplane 
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are used to designate airplanes flying, respectively, at speeds far below 
those at which dynamic and aeroelastic effects become important, and at 
speeds at which these effects have to be taken into account; the same 
airp lane can thus be "rigid" under some conditions and "flexible" under 
others . 

Several fundamental assumptions are inherent in the analysis con- 
tained in this paper. In the first place, all atmospheric disturbances, 
motions, and structural deformations are assumed to be small enough to 
produce forces that are linear and, hence, superposable. Also, the turbu- 
lent "input" to the airplane is assumed to be stationary in a statistical 
sense; that Is, the turbulence in the plane of the flight path is homo- 
geneous. For the large airplane, the additional assumption Is made that 
the turbulence is axisymmetrie with respect to vertical axes, a condition 
less severe than complete isotropy. The statistical characteristics of 
the turbulence are thus assumed to be invariant under a translation of 
the space origin within the horizontal plane and under a rotation of the 
coordinates about the vertical axis. Finally, Taylor's hypothesis to 
the effect that time displacements are equivalent to longitudinal space 
displacements Is assumed to be valid. 

The aerodynamic forces directly due to atmospheric turbulence, which 
constitute the input forces for the dynamic system represented by the 
airplane, are calculated in the first part of this paper for the large 
airplane, that is, for the case where the spanwise distribution of the 
Intensity of turbulence has to be taken into account. (The effect of 
spanwise variation of gust intensity on the lift has been treated by a 
slightly different method in ref. 5*) 

The dynamics of the rigid airplane are considered in the second part. 
The dynamic system Is now represented by a set of three simultaneous 
ordinary differential equations, rather than one aB in reference 2; none- 
theless, the problem of calculating the required transfer functions is 
still one of simple algebra. 

The next part is concerned with the small flexible airplane and thus 
has direct application to fighter-type airplanes and guided missiles 
operating at relatively high speeds. In addition to serving as a prelim- 
inary to the last part of the paper. The dynamic system Is now represented 
by a partial differential equation, and the calculation of the transfer 
functions requires the solution of ordinary differential equations. Once 
these functions are calculated, however, the statistical tec hni ques are 
the same as before, as a result of the fact that the lateral variation in 
gust intensity is ignored. Either modal or numerical- integration approaches 
may be used to analyze the dynamics of a swept-wing airplane with arbitrary 
stiffness and mass distribution. Although modal approaches have usually 
been preferred in the past for similar problems, it was believed that, in 
view of the highly complex nature of modern aircraft structures and the 
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advanced type of computing machinery required and generally available for 
their analysis, the numerical-integration approach would be preferable 
and it has, therefore, been used. 

The last part contains the analysis of the large flexible airplane. 
The statistical problem is now that of a system which is characterized 
by a partial differential equation with time and a space coordinate as 
independent variables and which is subjected to a random input that 
varies in time and space, so that more is required than the transfer 
functions from the gust intensity at one point on the wing to the stresses 
at another. The particular statistical problem presented by this case 
is considered in some detail, and the appropriate transfer functions are 
then obtained by using the numerical-integration approach presented in 
the preceding part to solve, in effect, the ordinary differential equa- 
tions which describe the wing deformations at any given frequency. 


SYMBOIS 


A 

b 

C(k) 

c L 

C l 

% 

Cmq 

Sx 

c 

c 


aspect ratio 
span 

Theodors en function 
lift coefficient, L/qS 

lift-curve slope 

rolling-moment coefficient, L’/cjSk 

coefficient of damping in roll, defined as positive for 
positive damping 

daznping-ln-pitch derivative 

static pitching derivative 

chord, parallel to plane of symmetry 
average chord, S/b 

section lift coefficient at station y, l / qc 
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El bending stiffness 

e-j_ distance from section aerodynamic center to shear center, 

fraction of chord 

eg distance from shear center to section center of gravity, 

fraction of chord 

e, distance from shear center to the midchord point, fraction 

^ of chord 

e> distance from shear center to the 3A-chord point, fraction 

of chord 

dimensionless lift-influence function (Green's function for 
the spanwise lift distribution) 

GJ torsional stiffness 

g acceleration due to gravity 

H(o>) response to sinusoidal oscillation, Fourier transform of h(t) 

h(t) indicial-response function 

mass moment of inertia about X-axis 

Iy mass moment of inertia about Y-axis 

Bessel functions of the first kind, order 0 and 1 

Kq,K]_ modified Bessel functions of the second kind, order 0 and 1 
Kig integral of Kq 

k reduced frequency, coc/2U 

k' dimensionless frequency, ojL*/(J 

L lift 

L 1 rolling moment 

L* integral scale of turbulence 

l distributed lift per unit distance along the span 
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M 

Mb 

% 

m 

a 


q. 

r 


S 

s 

t 

U 

u 

Up 

V 

W 

W 

V P 

X 

Ax 




pitching moment 
h ending moment 

twisting moment 

mass (of airplane., unless designated otherwise by subscripts) 

distributed mass per unit distance along the span 

distributed twisting moment (about axes perpendicular to the 
plane of symmetry; per unit distance along the span 

dynamic pressure 

radius of gyration about center of gravity; longitudinal dis- 
placement corresponding to time displacement r 

wing area 

T* 

scale parameter, 
time 

mean flying speed 

longitudinal component of gust yelocity 
horizontal component of disturbed motion 

lateral component of gust velocity 
weight of airplane 
vertical component of gust velocity 
vertical component of disturbed motion 

coordinate along mean flight path 

distance from intersection of elastic axis and root chord to 
airplane center of gravity 

tail length, distance frcm airplane center of gravity to 
aerodynamic center of tail 
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y 

z 


p 

r( tj) 
7(y) 
r'(y) 

n 

e 

A 

p 

T 

4 » 

* 

0(k) 

9 

9 

t 

¥ 


CD 


modified tail length, distance from. intersection of elastic 
axis and root chord to aerodynamic center of tail 

coordinate perpendicular to plane of symmetry 

coordinate in plane of symmetry perpendicular to mean flight 
path; vertical deflection 

inclination of chord to X-axis 

span ratio, b/L* 

autoconvolution function for /(y) 

dimensionless lift distribution, cc^jcC-^ 

dimensionless lift distribution in roll, cc^jcC-^ 

variable of integration corresponding to y 
angle of pitch 
angle of sveepback 
mass density of the air 

time displacement, argument of time-correlation function 
the power spectrum $ in the case of axisymmetzy 
two-dimensional power spectrum (double Fourier transform of If?) 
Sears function (unsteady- lift function for gust penetration) 
one-dimensional or point power spectrum 

two-dimensional power spectrum (single Fourier transform of ijf) 
one-dimensional or point correlation function 
two-dimensional correlation function 
frequency of oscillation 
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Subscripts : 


e,e 

effective 

f 

fuselage 

r 

wing root 

t 

tail 

u 

horizontal component of turbulence 

w 

wing; vertical component of turbulence 

Matrix 

notation: 

[ ] 

square or rectangular matrix 

n 

diagonal matrix 

L J 

row matrix 

{} 

column matrix 


Dots over symbols indicate derivatives with respect to time. 

AERODYNAMIC FORCES RESULTING DIRECTLY FROM 
ATMOSPHERIC TURBULENCE 


The motions of a rigid airplane depend on the overall forces and 
moments, whereas the stresses of a rigid airplane and the motions and 
stresses of a flexible airplane depend on the distribution of these 
forces, as well. This part of the report is concerned with the calcu- 
lation of the integrated and distributed forces and moments directly 
due to atmospheric turbulence when the spanwise variation of gust inten- 
sity has to be taken into account. (The forces and moments caused by 
the motions which result from the forces treated in this part can be 
calculated by conventional methods and will not be considered here, 
although the combined forces will be considered in the following partB.) 
Thus, this part serves as a basis for all the material presented in the 
later parts pertaining to the large airplane, and, hence, the fundamental 
notions required for an analysis of the large airplane are introduced 
here and discussed in seme detail. 

The basic approach is as follows: First, the instantaneous value 

of the quantity of interest, such as the lift, is expressed in terms of 
the instantaneous gust intensity at a point and a suitable influence 
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function. For the lift this influence function can and will he identified 
with a certain lift distribution on the given wing in reverse flow; the 
same procedure may and, in the case of the rolling moment, will he fol- 
lowed for other integrated forces. On the other hand, for the local lift 
the influence function is the Green's function for the three-dimensional 
unsteady- lift problem and cannot he identified with an easily calculated 
lift distribution on the wing in reverse flow. Inasmuch as no knowledge 
concerning this function appears to he available, a method of calculating 
an approximate Green's function for this problem is outlined herein. 

The required influence functions for integrated effects can be synthesized 
from this function, and if the associated lift distribution in reverse 
flow cannot be calculated . conveniently this approach may be preferable. 
This technique is illustrated here by means of the bending moment. 

The next step consists in using the expression for the instantaneous 
value of the given quantity to calculate a correlation function for this 
quantity in terms of a correlation function of the normal component of 
the atmospheric turbulence. The power spectrum for the given quantity 
can then be obtained by taking the Fourier transform of its correlation 
function. This power spectrum is considered herein to be the desired 
end result, because the mean-square values of the quantity and its deriv- 
atives can be obtained from it, and other statistical parameters of 
interest can be obtained from these mean-square values. Several alter- 
native approaches for calculating the aforementioned spectrum, either 
from the correlation function or directly from the spectrum of atmospheric 
turbulence, are given in connection with the lift and are directly appli- 
cable to other quantities as well. 

In this part of the report the assumption is made that the influence 
functions of concern can be written as products of a function of time 
alone and a function of distance along the span alone, and advantage is 
taken of this simplification in calculating the desired spectra. This 
restriction is removed in the last part of the report, where the inf luence 
functions considered cannot be separated into space-dependent and time- 
dependent constituents, and so the approaches outlined there are general- 
izations of those presented in this part; they may also be used for the 
calculation of the spectra of the quantities considered in this part if 
the assumption concerning the influence functions is not valid. 


Definitions of Statistical Parameters 

As pointed out in the introduction, the intensity of the vertical 
component of turbulence w(t) is a random process, so that the result ing 
forces can also be known only in a statistical sense. The purpose of 
this part is to calculate certain statistical properties of these forces 
namely, their mean-square values, their correlation functions, and their 
power spectra. The fundamental principles involved in statistical 
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analyses of the type considered herein are expounded in seme detail, and 
citations of the literature on the subject are given in references 1 
and 2. These fundamentals will therefore not be repeated here. However, 
both for the sake of ready reference and inasmuch as the statistical 
terms are not always defined in the same manner, the forms that are used 
herein are indicated in the succeeding paragraphs. 

The time average of a time -dependent quantity is designated by a bar 
placed over the symbol, 1 and is defined as follows: 



The assumption will always be made that this limit exists and is invariant 
under a translation of the origin of time. This assumption implies tha t 
the processes considered here are stationary in a statistical sense. 

The mean of a random process f(t) is defined as its time average, 
and is always assumed to be zero. In dealing with processes with nonzero 
mean this analysis is thus pertinent only to the process which consists 
of the difference between the original process and its mean value. Simi- 
larly, the mean-square value of a random process f(t) is defined as 
the time average of the. square of the process, so that 

f 2 = lim -i- / f 2 (t) dt 

T — > oo 2T y 

The time -correlation function of f(t) is defined as 


(t ) = f(t) f ( t+T ) 


so that 


= t f (o) 


when no possibility of confusion exists, a bar is also used to 
designate a space average, as in the case of c and y. Also, for the 
components of turbulence the correlation functions depend primarily on 
space displacements and can be defined by space averages. For the sake 
of consistency, however, they are considered to be defined by time 
averages . 
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the power spectrum of 
correlation function: 

f(t) 

is defined as 

the 

cp f (ai) 

II 

r e- toT *f(T) 

dr 


U 

—CO 



II 

-3S>r 

L * 

dr 


The second form is the more convenient one when i{f f represents more 

nearly a space correlation than a time correlation, so that it depends 
directly on a space displacement r = Ut rather than on a time displace- 
ment t. For the purposes of the analysis presented herein, the gust 
or input correlation functions have this property, so that, for instance, 
Tjr w may he defined as 


* w (r) = v(t) w(t+£) 


By virtue of the reciprocal properties of Fourier transforms and 
the symmetry of ijr f (r), which is a direct consequence of the assumed 

stationarity of f(t), f 2 may he expressed in terms of the spectrum of 
f as 



CD 


If f(t) is the Input of a linear system, the power spectrum of 
the response x(t) of the system is related to the power spectrum of 
f(t) hy 


<P x (oO = |H(aa)| 2 cp (co) 


( 2 ) 


where H(cn) is the transfer function of the system, that is, the complex 
amplitude of the response of the system- to unit sinusoidal input. Hence, 
H(co) is also the Fourier transform of the indicial response h(t) of 
the system, which is defined herein as the response of the system to a 
unit impulsive input: 



h(t) dt 


(3) 
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where the lower limit could "be taken as zero, since h(t) is zero for 
t < 0. Conversely, h(t) can he obtained from H(cc) by means of the 
inverse of equation (3): 

h(t) = — f E(oi) cks (3a-) 

2 * J _co 


The mean-square value of the response can then be obtained by inte- 
grating its spectrum. Similarly, the mean-square values of the nth time 
derivative of the response can be obtained from the (2n)th moment of the 
output spectrum. For instance, ■ - - 


(if = f 0 ' r > ) “ 2 

From the mean-square values of these derivatives other statistical quan- 
tities of interest, such as the expected number of peaks of the response 
per unit time, can then be calculated. 


Lift-Influence Functions in Unsteady Flow 

At any time the lift on a wing which results directly from atmos- 
pheric disturbances can be expressed for an unswept wing as 



where h(t,y) dy is a lift-influence function which represents the lift 
caused by an impulsive vertical gust of Width dy which at time t = 0 
impinges on the wing at station y. 

The Influence functions required in equation (4) are difficult to 
calculate directly} methods for obtaining lift distributions on wings of 
finite span in unsteady flow usually require numerical solutions which 
do not lend themselves readily to the analysis of angle-of -attack distri- 
butions represented by delta functions. However, by virtue of the reci- 
procity theorems of linearized lifting-surface theory (ref. 6, for 
instance) the lift influence function for a twisted wing in indicial 
motion is equal to the lift distribution on that wing during indicial 
motion in the reverse direction with uniform unit angle of attack. The 
lift distribution In indicial motion with uniform angle of attack can 
be calculated relatively easily. 
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For the few cases for which calculations have teen made (namely, 
same maswept wings), this lift distribution tends to be substantially 
invariant in time, except for overall magnitude. For instance, the 
calculations of reference 7 indicate that the lift distribution of an 
oscillating rectangular or elliptic wing in incompressible flow is sub- 
stantially independent of frequency, so that in indicial motion it is 
substantially independent of time. This simplification may not be valid 
for swept wings. 

The lift influence function can then be written as 

k(t,y) = i h L (t) r(y) (5) 

where y(y) defines the steady-state lift distribution for uniform unit 
angle of attack: 


r(y) = 



and where hr(t) describes the variation of the overall magnitude of 
the lift as a function of time after entry into a sharp-edge gust and 
may be written as 


h Ti (t) = 


C I*x' 


qS 


U 



In turn, kg is the lift response to a unit Bharp-edge gust normalized 
to a steady-state value of unity and is, as in the preceding equation, 
usually expressed in terms of semi chord lengths traveled The time 

derivative is taken here because the response wanted is the one to a 
unit impulsive gust rather than a sharp-edge gust. 

The Fourier transform E-^(cs) of this function h-^(t) is proportional 

to a function 0(k), which may be termed the generalized Sears function 
because for two-dimensional incompressible flow it is the Sears function: 

000 ( 6 ) 


where 


k = ^2. 

2U 
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The function H^cd) represents the complex amplitude of the lift due 

to sinusoidal gusts of unit amplitude and is thus the transfer function 
from the gust to the lift. 

Actually, only the absolute square of will be required. For 

two-dimensional incompressible flow the following approximation is given 
in reference 1: 


I0OOI 2 


1 

1 + 2iik 


(T) 


This expression has the advantage of simplicity, although it is somewhat 
in error compared with the abolute square of the Sears function at very 
low frequencies, a fact which could be remedied by using the approximation 


| 0«| 2 


1 + ah 

1 + ak(l + 2rtk) 


where a is about 15 for a good overall fit to the exact expression. 
However, the behavior of the Sears function itself at very high frequencies 
is unrealistic, because its absolute square goes to zero as l/k, whereas 
for any nonzero Mach number and any finite span the absolute square of 
the generalized Sears function can be shown to tend to zero at least as 

rapidly as l/k 2 . As a result of these discrepancies, the approximation 
given by equation (7) and, for the same reason, the absolute square of 
the exact Sears function, cannot be used to bbtain moments of the lift 
spectrum, that is, values of the mean-square derivatives of the lift, 
although they may be adequate for calculating the mean-square value of 
the lift itself in many cases, particularly when absolute accuracy is 

not required. Whenever possible the values ~of j 0(k) | given in refer- 
ence 8 for the plan form and Mach number of interest should be used. 


Mean-Square Lift and Its Spectral Resolution for 
the Unswept Wing 

Basic equations .- The correlation function of the lift can be 
expressed, by virtue of equation (4), as 


* l (t) J h ( t l' y l) h ( t 2' y 2) W ( U ( t ' t l) > y l) v(U(t*r-t2),y 2 ) ^ dy 2 it x dt 2 (8) 
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where the averaged product on the right side represents a velocity corre- 
lation function. This function depends in general on hoth space and time 
displacements. However , if Taylor's hypothesis is made, the time dis- 
placements are equivalent to longitudinal space displacements. The veloc- 
ity correlation functions are then functions only of longitudinal and 
lateral space displacements. Thus, for homogeneous turbulence, 

w(x+Ut,y) w(x+|+U(t+T),y+T)) = t)t w (4+[Jt,t]) ( 9 ) 

Tn addi tion to Taylor's hypothesis and the assumption of homogeneity, 
the turbulence is assumed to be axisynmetrlc with respect to vertical 

axes, so that ijr (g , tj) is a function only of p 2 + t] 2 . This function 

is the ordinary space-correlation function i|r , so that 

t v (UT+4,n) = * w (^(Ut+6 ) 2 +t} 2 ) 


and, hence. 


V t) = jl %r (l^ 2 ( T+t l" t fi) 2+ ( y 2" y l) 2 ) ^2 dt l dt £ 


If the assumption implicit in equation (5) is now made, the preceding 
equation can be written as 


/ CO p 00 

J h L{H) h L( t 2) v( u ( r + t l- t 2)) dt l dt 2 ( 1Q ) 

—CO ^ — OO ^ V / 


where 


W Dt) = r ( y l) 7 fr) ^ 


= r(t]) i w (ju 2 T 2 -t-r, 2 j (11) 

where, in turn, F(t]) is an autoconvolution of y(y) defined by 

Hi"" ^ 

rU) = |f / 7(7) /Cm) <3 y (12) 

* J -b /2 
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The validity of equation (ll) can he demonstrated hy performing the 
integration in the y-^yg plane as indicated in the following sketch: 



That is, integration is performed first over y-^, with a variable 
q = y 2 - y^ held constant. For this integration ^ is constant, so 
that only 7(y"i) and 7(y]_ +r l) are involved, and the result is a 

function of q which is one -half of the function r(q) defined hy 
equation (12). The second integration is then performed over tj, 
yielding equation (ll) except for a faetbr of l/2. In this process, 
only the part of the square above the line y 2 = y-j_ is covered. How- 
ever, by a similar process, the part of the _lntegral corresponding to 
the part of the plane below this line can be evaluated and shown to be 
equal to the first part, so that by defining T as in equation (l2), 
both parts are taken into account simultaneously in equation (ll). 

The quantity 



* Jj o r(n) t w (n) <3-n 

may be considered an averaged mean-square vertical component of turbu- 
lence; ijc^ (Ur ) is then the corresponding correlation function, and the 

Fourier transform of the latter. 
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v (a>) = is J e ^ e (ur) d ^ UT) (i3) 

is the corresponding power spectrum. 

Once cp^Cco) has been obtained, the power spectrum, of the lift 

can be obtained by taking the Fourier transforms of both sides of equa- 
tion (10). The result is 

<P L (oi) = |h l (co)| 9 Ve (<x>) (14) 

where Hj-^cd) is the transfer function defined in the preceding section. 

The mean square of the lift can then be obtained by integrating its 
spectrum, as indicated in equation (l). 

Equation (l4) has the same form as the corresponding equation for - 
the case where spanwise averaging of the effects of turbulence is not 
taken into account (see eq. ( 2 ) and ref. l), except that cp^fo)) is now 

replaced by cd (cjd). Thus, the spectrum of the averaged turbulence must 
w e 

approach that of the unaveraged turbulence when the span approaches zero, 
as may be seen to be the case from equation (ll) and the definitions of 
7 (y) and nCn) . 

Two alternative approaches . - The defining relations for <Pw e (cD), 

equations (ll) and (13), do not necessarily represent the best method 
of calculating it in any given case. A slightly different expression 
appears to be more convenient in general. It consists in substituting 
equation (ll) into equation (l3) and inverting the order of integration 
to yield 

c f t w e ( a) ) = b J 0 ?w^^) d, l U5) 

where ■ 

9 w (cd,T]) = ^ J e 4 (Ut) (l6) 

so that 

<P w (cu,0) = cp w (o3) 
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A third approach, which has certain advantages over the others, is 
similar to the one that has "been used in reference 5* I 11 this approach 
the assumption of axisymmetry is not made initially, and use is made of 
the spectrum 


®w( X 1^2) 



e _i ( A l^ +A 2 T l) 


<5.1 <3-n 


(17a) 


The correlation function ih. (Ut) 

w e 


can then be written as 


^ (Ut) = \ 


^ 7(y2) ^ ^ ®1 “2 


Substitution of this expression into equation (l3) and then interchanging 
the order of integration yields : 


= k f $ v(t^2)l ? ( x e)| 2 % 


where 


rh/ 2 

r(x) = i / r(y) ay 
b J-b/2 

and is, as a result of the symmetry of 7(y), real and symmetric in A. 

Now, if the turbulence is axisymmetric, 4^ depends only on 
A s ^A-^ 2 + A 2 2 > that is, 

®w(^1^2) = *(PY) 

so that the expression for <p w (oj) becomes 


(1Tb) 


V Ca)) = U 0 Mp) 2 ^ 2 )^ " 


(18) 



NACA TN 3910 


19 


where ® W (A) is related to i|r (r) "by 

r°° 

® W (A) = | J o r J cA r ) t w O) ^ 

Although the spectra used in the last two approaches hare been 

defined by expressions involving the point correlation function i|r (r), 

w 

they can be expressed equally well in terms of the point spectrum cp (cd), 

so that if, say, an experimentally obtained point spectrum is to be used, 
it need not be transformed into a correlation function before it can be 
used in these calculations. The required relations are 


%>,!}) = qvfa) -yf ^(“i) 

fo-L 2 - 




2 2 ’ 
CD-j_ —CD 


U 


(19) 


CD 


and 



( 20 ) 


where the notation 1 is used to specify that the finite part of 
the integral is to be taken,- an operation which may be performed by 
integrating by parts and ignoring the infinite part, so that, in terms 
of a proper integral. 




2U 2 P” ^wfc) dcp 

* Au fco ^2 . X %2 


Also, the function ® W (A) can be obtained frcm ^(cd, t] ) by means of 
either of the relations 


«„(*> 





= ^ COT ' n) ^ 
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Results of calculations .- In order to illustrate the magnitude of 
the effects under consideration , calculations have teen mad'fe for a uni- 
form. loading ?(y) = 1 and a point correlation function which has teen 
used in references 1 and 2 and appears to, fit experimental data (ref. 4) 
fairly well over a large portion of the significant frequency range, 
namely. 



so that 


<P w (a>) “ 


w 2 !* 

rtU 


1 + 3k* 2 
(l + k ' 2 ) 2 


where L* is the integral scale of turbulence, which is here defined as 


L* = 



dr 


and where 



This correlation function has the drawback that the moments of the 
spectrum associated with it are infinite, so that it Implies a process 
with infinite mean-square derivatives, but it is quite useful if only 
the mean-square value of the process itself is of interest. 

For uniform loading. 


and 


rW = 2(1 -£) 


- 


sin 


Ab 


Xb . 
2 
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and for the given correlation function, 

..2 




nU 


fe 1 + ^ - fef -^2 

* , w.2r' \ L / vl; L + k * 


I JL 

L 


f (i + k- 2 ) 

•where Kq and are modified Bessel functions of the second kind, and 


=■ 


. 2 T * 2 


-.2 


3 v^L* k 

■2 “ 372 


" (1 + k- 2 )' 

The mean-square average gust intensity for this case is given by 

2 ~2 1 - e“ P 

v = v 

6 3 

where p = b/L*, and is shown in figure l(a). The correlation function 
is given by 


W r) - f 


G ainh-1 |) - o 2 ^( g ; .1 nh-1 |) + d(e-“ - e-t 5 *?) 


where cr = r/L*, and where KqCctjQ) and K-j^crjP ) are incomplete modi- 
fied Bessel functions of the second kind defined by 


K^G) COBh ^ e ~* C ° Sh * 


d-a 


The spectrum qu (cd) is given by 


m / \ w 2 2L* 1 

^ p2 JtU 2N 3 

P (l + k* 2 ) 


3k ,2 p^Jl + k' 2 |kl 0 (p ^L+k’ 2 ) - 


p||l + k* 2 KQ^Jl+k’ 2 ) + (l - 3k’ 2 ) [2 - 2 pjl + k’ 2 K 1 (p^l+k’ 2 ) 

p 2 (l + k’ 2 ) Ko^+k’ 2 ) 
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where K1 q(x) is the integral of K^(x): 

^(xi) dx-L 

The functions ^ and cp v (normalized, with the averaged mean- 

square turbulent velocity) are shown in figures 2 and 3, respectively. 

The effect of the span ratio on the normalized correlation function of 
the averaged turbulence may be seen to be relatively small, so that the 
effect on the unnormalized correlation function is primarily the decrease 

in overall level given by the ratio of w~^ - to w^. A similar statement 

may be made for the power spectrum. If the power spectrum were not 
normalized the averaging effect of the span would tend to reduce the 
intensity of the spectrum at all frequencies, but the high frequencies 
would be attenuated much more than the low ones, as might be expected. 

In fact, although the unaveraged spectrum decreases as co at high 

frequencies, the averaged spectrum decreases as . 

The asymptotic values for b/l* — are shown in figures 2 and 3 
in order to indicate the nature of the functions considered here when 
the scale of turbulence is small compared with the span, as may be the 
case for a wind-tunnel model responding to natural or artificial tunnel 
turbulence, or for a buffeting wing or tail surface, although this con- 
dition is not of practical concern for the gust-load problem. 

The power spectrum of the lift is equal to the product of the power 
spectrum cp w ^(ca) and the absolute square of the transfer function H^(oi) , 

as indicated In equation (l4). Inasmuch as this lift is not an end in 
itself but only one of the parameters that enter Into the calculations 
of the motion of the airplane, its mean-square intensity Is of little 
practical significance} its spectrum is the quantity needed in further 
calculations. However, If the mean-square intensity is wanted for any 
reason It can be obtained by integrating the spectrum. Thus, for instance, 

the approximate expression for |0(k)| given by equation ( 7 ) and the 
spectrum 9 w (o>) used for the preceding calculation yields the mean-square 

lift: 
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where 


s 



and, in view of the observation that much of the turbulent energy is 
contained in a region for which the span has a very small effect on the 
(normalized) spectrum, this equation should serve to furnish an approxi- 
mation to the mean-square lift for nonvanishing span, provided is 

used instead of w^. However, the mean-square values of the derivatives 
of the lift cannot be calculated in this simple manner, because the 
deviation of the normalized averaged spectrum from the point spectrum at 
high frequencies cannot be ignored in calculating the moments of the 
lift spectrum. 


Mean-Square Lift of the Swept Wing 

For the yawed or sideslipping unswept wing, equations ( 10 ) and (l4) 
for the lift-correlation function and spectrum are still valid if an 
appropriate lift-influence function is used, and if the correlation 
function for the averaged turbulence is defined by 

pb cos A / p — — \ 

* We (UT) * b "cos ' A J 0 ^w\^ U t+t 1 t£m A ) +T rJ ^ 

where r( iq) now pertains to a lift -distribution function 7 (y)* which is 

appropriate for yawed motion and is defined for - cos A = y = jj- cos A. 

The mean-square averaged intensity of the vertical component of turbulence 
is then 


— - _ pb cos A , _ . 

w e 2 - f r(T]) %(—~) dT) 

b cos A Jq ™ \cos A/ 


r(n'cos A) ct-T)' 


where ti * = — 3-—. Thus,. this mean-square intensity is unaffected by the 
cos A 

yawing process, except for 'the slight change which results from the change 
in r(ti), although the spectral resolution of the averaged intensity 
changes in the process. 
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For the swept. . wing both y 2 - yq and | y 2 | - |yq| occur in the 
integral, so that the reduction of the 'double integral for ‘ilf Tr (Ur) to 

w e 

a single integral (see eq. (ll)) cannot be affected so simply. The 
double integral for the swept wing is 

W UT > = & Xb/ 2 fJ/2 7 ( y 0 r ( y 2) ^[ UT+ (l y 2li y i| ) taa A ] 2 + (y 2 -yi) 2 ) ayi dy 2 (21) 

From this integral, by using rectangular lift distributions and the 
aforementioned point correlation function, the mean square of the averaged 
turbulence as well as the corresponding correlation function and power 
spectrum have been calculated by numerical integration for various sweep 

n 

angles A, the ratio — - — being maintained at 0-5- (The decision to 

cos A 

hold — ^ — rather than 8 constant was reached as a result of the 
cos A 

foregoing analysis, of the yawed unswept wing, which indicated that the 
effects of sweep should be minimized in this manner.) The results for 

w e 2 are shown in figure 1, and the effect of sweep on w e 2 iB seen to 
be small for this comparison. The calculated correlation functions and 

spectra (normalized with respect to w e 2 ) are not shown because they 
agreed with those for A - 0 within less: than 1 percent for most values 
of Ur and h’j respectively. 


Mean-Square Rolling Moment 

In the preceding sections the averaging effect of the span has been 
shown to consist, essentially, in reducing the effective intensity of the 
turbulence sensed by the wing; thus, it only modifies the forces present 
on a wing of small. span. If the analysis is" extended to the rolling 
moment, however, a new phenomenon appears". When a wing is so small relative 
to the scale of turbulence . that at any instant all of its points experience 
the same turbulent velocity, the wing experiences no rolling moment as the 
result of the direct action of turbulence, (although it may experience a 
rolling moment indirectly as a result of the rolling and yawing" motion 
caused by the lateral component- of the turbulence) . On the other hand, 
on a large wing the different intensities, of the turbulence at different 
points on the span give rise directly to a net rolling moment, which then 
results in rolling motion. In this section the. mean-square value of this 
moment is calculated. 

At any instant t the rolling moment L'(t) can be written in the 
same form as the lift L(t) in equation (!{•); however, according to the 
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previously mentioned reciprocity theorem, the lift-influence func- 
tion h(t,y) is now the lift distribution for an indicial roll with 
unit he lix angle at the wing tip. If the assumption of invariance of 
this distribution with time is made, as for the symmetric case (see 
eq. (5)), then the required lift-influence function can be written as 

h(t,y) = i h’(t) 7 ’(y) 


where 

7 '(y) 


Cu (t)qSb 

h‘(t) = — 2 - — ■ , and where the steady-state lift distribution 

3 J 

cc 7 

= = —^ now pertains to a unit linear antisymmetric angle of attack. 
cC 1 

The correlation function for the moment can then be written as 


a 03 

h '( t i) “’(t?) dt i at 2 < 22 > 

00 


where 


^e' (aT) 7 '( 3r l) f)^l ^ 


- J J r ' Cti) t w (ju 2 T 2 +T] 2 ) dTj 


(23) 


where, in turn, in analogy with equation ( 12 ), 

o( b/ 2 ) -q 

r'U) = f / , 7'(y) 7'Cy-Hi) ay 

b vi-b/2 

Hence, the second and third approaches indicated in the section concerned 
with the mean-square lift (see eqs. (15) and (16)) can be used to obtain 
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and 


where 


and 



r’(A)| 2 dA 


?'(A) 



“ (“Sr") l ^ k ^| 2 


(25) 


so that the mean-square rolling moment can. he obtained by integrating 
this spectrum. 


A qualitative indication of the effect of span on the mean-square 
rolling moment sensed by an airplane may be” obtained from the quantity 

ik_ r(o), which represents the integral of the spectrum cp *. For a 
w 6 w e 

linear loading, y' - 6 and 


r'Crj) 



r'(A) 


6i 


Ab 

IT 


cos 



Ab 

2 



Hence, using the aforementioned expression for the point correlation 
function yields 


/ \ _ i/v 2 18 

9 W « Coo; = 


itU . n .k 

P'^fl + k' 2 ) 


[-32 + 2P' 2 + (l6p’ 2 + 3' 4 ) KqCp') 


( 323 ' + 6p' 3 ) k l ( 3 ') + k’^ J 32 - 6 p ' 2 - 163' 2 Kq( 3 1 ) - 


( 323 ' + 23 ,5 ) k l ( 3 ') 
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where 


p' = p](l + k’ 2 

and 

* We ' (°) = " 6 + e" P ( 6 + 6p + 3|3 2 + £ P 5 ) 

This spectrum and its integral approach zero linearly as p tends to 
zero. 


Generalized Aerodynamic Influence Functions in Unsteady Flow 

The aerodynamic influence functions used in the preceding sections 
define the contribution of a given station of a wing to the total lift 
and rolling moment. In the analysis of a flexible wing, and even in the 
calculation of certain properties of a rigid wing, generalized aerodynamic 
influence functions are required, which define the contribution of one 
station on the wing to the lift at another station and thus represent a 
Green's function for the unsteady spanwise lift distribution. No work, 
appears to have been done on such functions. For steady flow, apart from 
seme calculations for supersonic speeds which are based on the subdivision 
of a given wing into a number of squares, the only available results 
appear to be those given in references 9 aud 10. 

The analysis in this section is based on reference 8 and consists 
in a generalization of the method presented therein to -unsteady flow. 

This method constitutes an attempt to predict the lift distribution for 
any given twist on the basis of knowledge concerning a few definite 
angle-of-attack distributions, and may therefore be termed a function- 
interpolation method. For the present purpose, the presumably known lift 
distributions are the ones for uniform angle of attack in direct and 
reverse flow in a dimensionless form, namely, cej/cC^j they will be 

referred to as 7jj(y) and 7 j^(y}, respectively. (The function 7 (y) 
used previously is the one now designated by 7p(y)») Also required Is 
the coefficient of damping in roll (The lift-curve slope and the 

coefficient of damping in roll are the- same in direct and reverse flow 
by virtue of the reciprocity theorem, so that no distinction will be 
made. ) 

The approach of reference 9 then yields the following approximate 
expression for the lift distribution due to any angle-of-attack distribu- 
tion: 
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-=^ = Cj^j a 7 D (y) + x[a(y) - a] y R (y)| 


(26) 


where 


r b/2 

a = »J- b/2 7E(y) “ (y i 


<3y 


°Zp 
K = — £ 


b5 


C Ia rb/2 


(27) 


p D/ 

4 / 7 R (y) y 2 ^ 

■j 0 


Values of K may be obtained from the information given in references 9 
and 10. As the aspect ratio tends to zero, K approaches l/2, whereas 
for aspect ratios approaching infinity, K tends to 1. The following 
relations can be obtained from elementary definitions and from the 
aforementioned reciprocity theorem: 




With the aid of these relations the lift distributions given by equa- 
tion (26) may readily be seen to have the correct lift and rolling moment 
for all angle -of -attack distributions, and to reduce to the exact lift 
distributions for angle-of-attack distributions which vary linearly along 
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the span. A lift distribution which possesses these properties could 
readily be obtained by approximating any given angle -of -at tack distribu- 
tion by a suitably chosen linear one. However , this procedure would 
yield poorer approximations than the ones furnished by equation (26); 
for a parabolic angle-of -attack distribution on a wing of very small 
aspect ratio, for instance, this procedure would yield the lift 
distribution 




with a bending moment too low by 17 percent compared with the one of the 
exact lift distribution. 


^ ■ If ' tisf 


whereas equation (26) yields 

CC, A 

-T-tf- 


{&) 


1 + 2 {- 3 L 
,W 2 J 


1 + 


with a bending moment 8 percent too high. For very large aspect ratios, 
7 D and both tend to the chord distribution c/c, so that the lift 

distributions given by equation (26) tend to the correct limiting value. 


cc 


_L 

c /A- 


= C, g a 
a c 


(28) 


Therefore the accuracy of the results furnished by equation (26) may be 
expected to increase as the aspect ratio increases, whereas the accuracy 
of the other approximation is independent of aspect ratio. 

In this connection, it may be mentioned that "strip theory" consists 
in using equation (28) for all aspect ratios and, hence, is not very 
satisfactory for wings with medium and low aspect ratios. For instance, 
for the case discussed in the preceding paragraph, it furnishes a bending 
mcment which is too low by 25 percent for a delta wing, and too high by 
25 to 100 percent for tapered wings. (The two figures pertain to taper 
ratios of l/h and 1, respectively. ) 

Substituting the value for ' a from equation (27) into equation (26) 
gives the following expression for cc^/c: 

-=r = f jjjbCy) " K 7 R( y ) + bK s (y- T i)J 7 R U) j °,(ti) 
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The expression in braces in the integrand is the desired generalized 
aerodynamic influence function and -will be designated by G(y, tj),' so that 

o(y,T]) = [/ D (y) - K 7 R (y) + bK 5(y-n)J 7 r (t)) (29) 


and 



( 50 ) 


In the limiting cases of wings of very low or very high aspect ratio, 
and 7 r approach a common value, say 7. Thus, for wings of very 

low aspect ratio. 


o(y,n) = ^[/(y) + 13 s(y-t))J rU) 


( 51 a) 


where 6 is the Dirac delta (unit-impulse) function, and for wings of 
very high aspect ratio, 

G(y,r]) = b 6(y-i}) 7(1)) (31b) 

which is the Green's function associated with strip theory. Thus, on 
wings of very high aspect ratio all the lift produced by the local angle 
of attack at a given station is carried in the immediate vicinity of 
that station, whereas on wings of very low aspect ratio much of the lift 
Is carried elsewhere. This tendency for a given station to affect a 
greater portion of the wing as the span decreases is, of course, to be 
expected. 

The preceding analysis can be applied to the oscillatory case at 
a given reduced frequency k, as well as to the steady case. If the 
assumption of invariance of normalized distributions with time or fre- 
quency is made, as before, then 7p(y), 7 R (y), and K are independent 

of frequency, so that equation (30 ) can be written as 

Cr(k) qS p t/2 

*(y>k) - 5 / , O(y^) o,(tj) drj ( 32 a) 

b J -b/2 

where G(y, t)), defined as before, is independent of k. By applying a 
Fourier transformation to both sides of this equation the following 
relation is then obtained for flight through continuously varying tur- 
bulence (cf. eq. ( 4 )), as modified by the' assumption stated in equa- 
tion (5): 
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>b/2 


p n oi d 

dt l ^ G(y,Tj) w^Uffc-tJjTi) dT] (32b) 


where i(y,t) is the lift per unit span at station y and time t, and 
where the function h^(t) is the one used previously. 


The correlation function for this lift can then be written as 


W t) “ ^££m*i) MV) dt i at 2 $f' h/ lf* h/ l V’\) °( y ^) ^fVv^Vh) 2 ) d h dT >2 
~ \ [ f . h L Ch) ^(*2) dt l dt 2 ^ € ( U ( T+t l"V) ,y ) 

"h » — ca ^ —co 


(33a) 


where the function 

+v £ ( Ut f °(y^i) G (y^2) ^w(| u2 ' r2+ ( T1 2- T ii) 2 ) dT1 l d ^2 

represents an effective correlation function, which when transformed into 
the equivalent power spectrum q> w (ao,y) can be used to obtain the power 

spectrum for i(y, t) and hence its mean-square value. Thus 

<Pi( y )(®) |0(k)|\ e (^,y) (33b) 

Before the calculation of cp T ^(co,y) is discussed, the function 
*w e (tlT;y) will be defined in a somewhat more general form than in the 
preceding paragraph in order to anticipate future needs, namely 

V £ ( Ut ^1^2) f ^f 2 f 2 G ( y l ,Tl l) G ( y 2^2) %r(i u2T2+ (^l) 2 ) d \ d7 l 2 

If the assumption is now made that 7 ^(y) 7p(y) 8X6 H e same, as 

is the case for all unswept wings and for wings of very low or very high 
aspect ratio, they can both be identified with the function j(y) used 
previously, so that 

G(y,T)) * {(l - K) 7 (y) + bK 5(y-T))] /(tj) (3*0 
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( Ut ^i^ 2 ) = j ^ 1 “ Kj^CUr) + K(1 ' K) [^ e *( UT ' y l) + 

V*^ UT,y 2 )] + K \(fo 2 r H y 2 - y l f )} V 7 ( y l) 7 ( y 2 ) 

where iV„ (Ut) is the correlation function calculated previously for 
w e 

the averaged vertical component of turbulence, and where 

f We *(Ur,y) - if 7 ^ ^(iJu^+Cy-Ti) 2 ) df) 

Similar ly, then, ■ . ' 

^(“^1^2) = + K(1 “ K) ]%*(“ ,y l) + 

V*(“' y 2)] + k 2 9 w (^,y 2 -y L )| r(y L ) 7 { y 2 ) (33) 

where cp w ^(a),y) is the Fourier transform of t-vr e *(UT,y), so that 

r b/2 

cp *(oj,y) »i / /(t]) ?„-(^|y-T|) dl l (56) 

e D J--b /2 

and thus represents an averaged form of the spectrum cp^ introduced 

previously. (See eqs. (l 6 ) and (19).) For uniform spanwise loading and 
the point correlation function used previously, this function is 




L^w 2 1 

*U p 


(l + k' 2 )' 


1 ~5^fk , 2 [Ki 0 (a 1 ) + Ki 0 (a 2 ) - 


ai Ko( a i) - &2 M 8 ^)] + [Al 8 !) + j 


where 


— ± y f 

a i ' a 2 ^ Y + k ’ 


and is shown in figure 4 for several values of 
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The more restricted form of cp^ required in equation (33b) can 

now he obtained from the more general form given in equation (35) by- 
setting y 2 = y^ = y in the latter, to yield 




- (1 - 


K) 2 cp, 


w, 


>) + 2K(l - K) cp^*(co,y) 



(3T) 


Mean-Square Bending and Pitching Moments 

When the variation of the gust intensity along the span is taken 

- 5- 

into account, the mean-square lift I r is not an adequate index of the 
stresses in the wing, nor can the stresses be obtained from the mean- 

square lift distribution 2 2 (y ) . Instead, the mean-square bending and 
twisting moments, as well as the mean-square vertical shear, must be 
calculated directly. 

In a manner analogous to that employed for the lift and rolling 
moment, each of these quantities can be expressed in terms of a certain 
influence function which, by means of the reciprocity theorem, can be 
related to a lift distribution on the wing in reverse flow. For instance 
for the root bending moment the desired lift distribution is the one for 
an angle of attack which is zero on one wing and proportional to y on 
the other, as may be seen from the fact that the bending moment can be 
expressed as 

“b ’/-I,/ 2 f(y) ^ 

where 

f(y) = 0 (y < 0 ) 

f(y) » y (y ^ 0 ) 

so that, according to the reciprocity theory for un swept wings in steady 
or indicial flow. 



ct(y) <3y 


Similarly, for the root shear the required lift distribution corresponds 
to an angle of attack which is zero on one wing and uniform on the other. 
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Such lift distributions can be calculated readily. For instance, 
for subsonic flow and unswept wings they can be obtained from those given 
in reference lOj the lift distribution for the root bending moment is 
one-half the sum of a linear symmetric and linear antisymmetric lift 
distribution, and the lift distribution for the root shear is one-half 
the sum of a lift distribution due to a uniform angle of attack and a 
1 if f, distribution due to deflection of a full-span aileron. 


However, in seme cases such calculations may be time-consuming, 
anri an alternative approach may be desirable. One such approach consists 
in synthesizing the desired influence function from the generalized lift- 
influence (Green's) function discussed in the preceding section. In 
order to illustrate the use of this approach, it is adopted in this 
section. 

The bending moment at any station y ^0 ^ y ^ and at any 

time t can be obtained from the lift distribution considered in the 
preceding section as 


r b/2 

(y/t) = / (y 1 - y) i( y',t) dy 1 

J y 

pb /2 o 00 _pt/2 

- (y 1 - y) <3y' hj^t-L) dt x ij G(y',Ti) w^t-t^dn 


Mi 


B 


>b/2 


= b J dt l J h j 2 w( U (t-t 1 ),T|) dT) 


(38) 


where the influence function for the bending moment is 


>b/2 


^ p V / CL r- . 

M(y,n) = £ J ( y ’ - y) G(y’,Tl) dy' 


(39a) 


so that, upon Introducing the previously used function for G(y, t]) 

M(y,ri) = [(1 - K) M L (y) + K M 2 (y,Ti)] 7 (T 1 ) (53b) 



NACA TN 3910 


35 


where 


M-^y) 



- y) 7 (y') < ay' 


M 2 (y,n) = n - y (r» > y) 
MgCy^fi) =0 (*n = y) 


and, hence , for the root bending moment, 

vb/2 

M- 


1(0) = §J o y* r(y’) ay 1 

M 2 (o,n) = n (n > 0) 

1^(0, n) =0 (ti < 0) 


Hence, the correlation function for the root bending moment is 

p“ p“ nb/2 pb/2 f 

W T) MM MM 4t ^ 2 £J_ b/2 J. i/2 { (1 - K) V^ 

K(1 - K) M 1 (0)[M 2 (0,T ll ) + M 2 (0,T! 2 j] + 

K^(0,m) M 2( 0 ^)}y( T ll) rffe) v(j u2 (T+t 1 -t 2 ) 2 + (T12-T1!) 2 j dT] L dTi 2 

(40) 

The approach used in the preceding sections can now be used to obtain the 
power spectrum of the root bending moment by evaluating the inner pair' of 
integrals of equation (40), taking the Fourier transform of the result 
with respect to r, and multiplying the power spectrum obtained in this 

manner by — ) 1 0(k)[ 2 

indicated for the lift. 


or by using one of the alternative approaches 
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The spectrum corresponding to the inner pair of integrals is, for 
K = 1 and -uniform loading. 


Cp w . h(oj) = 


L*v^ 


rtU0 


(l + k' 2 ) 


(_6l4.-t. 20 12 ) + 80 ' 2 K o (^-) + (320’+ 0’ 5 ) K l ^) 


k' 2 |p ,3 ki o(y) + (64 - 60 ,2 ) - 80 ,2 K 0 (^) - (320 ' - 0 ,3 ) 


r 2 " 

= 0yl + k' , as before. 


where 0 

The integral of this spectrum is 


ilv -(°) = ^ 
30 


(-24 + + e“ P//2 (24 + 120 + 30 2 ) 


Examination of n(0) indicates that, although the mean-square bending 

moment tends to decrease as the span increases, it decreases less rapidly 

than the mean-square lift, with the result that the effective lateral 

center of pressure moves outboard. Quantitatively, the distance from the 

plane of symmetry to the effective lateral center of pressure can be 

defined as the square root of the ratio of the mean-square root bending 

moment to the mean square of the lift on one wing, that is, of the root 

shear. Although these mean squares have not been calculated, the square 

root of the ratio of iL, n(0) to the corresponding value for the root 

w e 

shear increases by 15.5 percent as 0 increases from O to infinity, with 
much of the increase realized at fairly small values of 0. 

For a swept wing the variation of the gust intensity along the span 
results in a pitching moment which must be taken into account in calcu- 
lations of the dynamic response of the airplane to continuous turbulence. 
This pitching moment can be obtained in substantially the same manner 
as the bending moment. Thus, if y is the station of the mean aero- 
dynamic chord, 

r b/2 

M(t) = tan A / (y - |y|) l(y,t) dy 
J -b/2 


= tan A 


_ p w p b/2 

b J MM dt l J M 3^) w ( U ( t- M' n ) <3.r) 


-b/2 


(41) 



NACA TH 3910 


37 


where 


ph/2 

M 3 (ti) o I / (y - |yf) G(y,Ti) dy 

J b J -b/2 


so that, with the previously used approximation to the Green's function, 
M 3 (ti) = |(1 - K)[y - 2M X (0)] + K(y - h|)| 7 (t0 


The correlation function, spectrum, and mean-square value of the pitching 
moment can then he obtained in the manner used in the preceding sections. 


Wing-Tail Correlation Effects 

The tail strikes a given gust some time after the wing does; as a 
result, a pitching moment arises which does not exist in steady or quasi- 
steady motion, nor if the airplane is very small, because then the time 
lag is insignificant. This pitching moment can, for the purpose of the 
present paper, be analyzed either in terms of the correlation between the 
gusts at the wing and those at the tail or, if a time-lag term is included 
in the indicial-response function of the tail, in terms of the correlation 
between the wing and tail response functions. The first point of view 
serves to exhibit the effect under consideration more clearly and is 
adopted first; the second is more convenient and is adopted in the subse- 
quent parts of this report. 

In the somewhat artificial case of a small wing and tail separated by 
a relatively large distance, only the distribution of turbulence along a 
line (the flight path) rather than In a portion of a plane is needed. 

The pitching moment due to the vertical component of atmospheric turbulence 
can then be written as 


M(t) h^t-J v(u(t-ti)) dt L + J h^t-J w(xf«J(t-t L )) dt L 


where is the tail length, and where h^t) and h.j.(t ) are the 

pitching-moment responses to judicial gusts hitting the wing and tail, 
respectively, at t = 0; both may Include unsteady-lift effects, and, 
if downwash effects are to be considered, h vr (t) should include the 

contribution to the pitching moment of the tail lift caused by the down- 
wash at the tail associated with the lift on the wing which results from 
the indicial gust. 
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The spectrum of this moment can then he written as 

q^(co) = | B^Ccd) +■ H t(cu)| 2 - 2R - (l - H t (aj) ' 9 v (o>) 

where the symbols R j | and * designate., respectively, the real part 
and the complex conjugate of a complex number. Obviously, when x^. 

approaches 0 the second term in the bracket in this equation vanishes, 
so that the first term represents the perfect-correlation effect, and 
the second represents the correction for Imperfect correlation. 

In order to furnish an estimate of the magnitude of the effects 
under consideration, same calculations have been made on the basis of 
the assumption that a real coefficient p exists such that 

H^ai) = p Hj-Cco) 

which implies that the attenuation with frequency of the contributions 
of the wing lift and tail lift, respectively, to the pitching moment is 
the same. The ratio p is -1 for neutral stability, and p > -1 for 
stable flight; it is positive when the aerodynamic-center location 
(tail off) is behind the center of gravity, so that positive values of 
p Eire not likely to be incurred with normal configurations and flight 
conditions . 

For this case, 

%(oi) = KM 2 [(1 + H) 2 - Zp(l - cos Sjj cp w (cu) 


tail alone 


= (1 + p) 2 - 2p ^ 


r*t *t 


|H t (cs)| 2 cp v (co) ^1 - cos da 

p °0 g 

J |Hfc(®)| cp w (cs) dm 
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The function f has been calculated for several values of its arguments 
by means of the lift -attenuation function given in equation ( 7 ) and the 
point spectrum used in the preceding sections, and is shown in figure 5* 
Also shown is the ratio of the mean-square moments as a function of the 
factor (J. which in figure 5 is referred to as The effect of 

imperfect correlation is seen to be very large as the condition of neutral 
stability is approached; the entire pitching moment is then the result 
os instantaneous differences in gust intensities at the wing and tail. 


In general, however, the tail length and the span are of the same 
order of magnitude, so that an analysis of the effect of imperfect corre- 
lation between the wing and tail must take into account the averaging 
effect of the wing span. The pitching moment at any instant Is then 


,b/2 


T P P °f 

M(t) = ^ / Wtq) dt i / y(y) w(u(t-ti),y] dy + 

^ -co -b/2 

r 00 

J ^t^l) ‘ W ’(x t + U (t-t x ) , 0^ dt L (42a) 


Hence, the power spectrum of this moment is 

<P M (c») = [Hy-Ccu)! 2 qv e («>) + | Kt-Cco)! ^(oj) + 

aR-Je 1 ^ Ht(o>) qv e *(«>, 0 ) ( 42 b) 

where <P We (co) is the averaged spectrum of equations (ll) and (13), and 
cp We *(cD,0) is the spectrum of equation ( 36 ) for y = 0. (it should be 
noted that h^(t) is the response to an indicial response which strikes 
the tail at t = 0; if it were the indicial response to a gust which 

i&ct 

strikes the wing at t = 0, the factor e u in the preceding equation 
would not be required. ) 


DYNAMICS OF THE RIGID AIRPLANE 


In this part of the report the motions of a rigid airplane subjected 
to atmospheric disturbances are considered. This part thus applies 
directly to those airplanes which fly at relatively low speeds and do 
not experience any significant structural deformations, and also serves 
as a preliminary to the treatment of the flexible airplane in the later 
parts . 
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The assumption is made that the motions are small enough to permit 
the use of linear approximations to the resulting aerodynamic forces and 
to permit the linear superposition of these ? orces . The longitudinal 
degrees of freedom (pitching * vertical* and horizontal motion) and the 
lateral degrees of freedcm (yaw* sideslip* and roll) can therefore he 
considered separately. The first two sections of this part are concerned 
with the longitudinal motions of a small airplane. The material presented 
here is thus a generalization of the single -degree -of -freedom analysis 
presented in reference 2. Mathematically* the dynamic system is now 
described by a set of simultaneous ordinary differential equations rather 
than a single one, but the problem of calculating the pertinent transfer 
functions is still one of simple algebra. The extension of these results 
to the large airplane is effected in the third section* using the tech- 
niques developed in the first part of this report. The lateral motions 
of a large airplane are considered briefly in the fourth section* and the 
combination of the stresses due to longitudinal and lateral motion of a 
large airplane is discussed in the fifth section. 


Equations of Longitudinal Motion 

The equations of motion of an airplane can be expressed in several 
coordinate systems . The system which is generally the most convenient 
one for analyses of airplane stability and is generally referred to as 
stability axes, consists of body-centered axes which are normal and 
parallel to the relative air velocity and rotate with it as the airplane 
pitches or yaws. (See ref. 11* for instance.) The aerodynamic forces 
related to this axis system can be measured more readily in wind tunnels 
than those related to other axis systems. In view of the very close 
relation of a stability analysis to the problem considered here* these 
axes will be used in this part of the paper* but in the analysis of the 
flexible airplane in the subsequent parts of the report space-centered 
axes will be used* because they are slightly more convenient for that 
purpose. 

The airplane will be considered to be in steady level flight prior 
to disturbance. The motions studied will be the deviations from their 
mean values; for instance* the angle 0 considered here will be the 
difference between the disturbed and the initial value of the angle of 
pitch. Hence* the motions and forces calculated by the method indicated 
here must be added to their mean values to obtain the total motions and 
forces . 

Inasmuch as* for the purpoBe.of a statistical analysis* the dynamic 
characteristics of the airplane are represented most conveniently by its 
transfer functions* attention will be confined in this section to sinus- 
oidal gusts and motions. 
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For this case the linearized equations of longitudinal motion can 
be written as follows (see eqs. 11-193 of ref. 11, for instance): 



-Zu 

-io)U 


Wp 


Zw 

Zu 

“Xw 

in) - X u 

g 

< 

Up 

> = 0(k) 

Xy 

Xu 

-ICEMjJ. - My 

-Mu 

-a) 2 - iuMq 


, 9 . 


Mw 

Mu 



(43) 


The stability derivatives which appear in these equations are defined 
in table 1 in terms of conventional aerodynamic coefficients, and the 
numerical values are given for the example used in reference 11. (The 
value of the mean chord is not given in ref. 11 but it is assumed herein 
to be 10 feet on the basis of other information given in ref. 11.) 

In analyses of the stability of a rigid airplane the quasi-steady 
approximation to unsteady-lift effects is -usually made, in which, in 
effect, the forces corresponding to a steady attitude, to constant dis- 
turbance velocities, and to constant accelerations are considered. This 
approximation is justified because the motions of concern are generally 
sufficiently slow. For the same reason this approximation can also be 
made in analyzing the response of an airplane to atmospheric turbulence. 

However, in this problem another type of unsteady-lift effect occurs, 
namely, that related to the forces directly attributable to the turbu- 
lence. This effect is here taken into account by multiplying the quasi- 
steady values of the forces due to gusts on the right side of equa- 
tion (43) by the attenuation function 0(k). This procedure implies the 
assumption that the airplane is small relative to the scale of turbulence, 
inasmuch as no averaging effects have been taken into account; these 
effects will be discussed presently. Also, this attenuation function is 
strictly applicable only to the normal forces. 

The unsteady effects on the drag are not known because of the rela- 
tively complicated nature of the mechanism which gives rise to drag. If, 
however, the assumption is made that upon entry into a sharp-edge gust 
the drag rises linearly and attains its steady-state value in the time 

required to travel N chord lengths, the drag equivalent of ] 0(k)| 2 
is the function 


1 - cos 2Nk 
2H 2 k 2 

which, for N equal to about 5 or 6, agrees fairly well with j0(k)J 2 
in the region of main interest (k > l). 
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The unsteady moment is also difficult to predict because of the 
paucity of knowledge concerning unsteady dcwnwash effects for wings of 
finite span. However, inasmuch as the wing lift contributes part of 
the moment and, through the mechanism of dcwnwash, determines to a large 
extent the moment contributed by the tail, the use of the lift attenu- 
ation function for the moment appears reasonable for a first approxima- 
tion, and the use of the same function for the lift, drag, and moment 
facilitates the analysis. 

For a more refined analysis, the time lag between the instants at 
which the gust hits the wing and the taii must be taken into account, 
not only for large airplanes, but even for small airplanes if phugoid 
motions are important. (See ref. 12.) Also, the lag in dcwnwash should 
be taken into account for large airplanes and possibly also for small 
airplanes in some cases . One way of achieving this result is indicated 
in the next part of the present paper; another is discussed in refer- 
ence 12. 

In equation (43) the unknown quantities Wp and Up are the normal 
and axial components of the disturbance velocities of the airplane rela- 
tive to the free stream. Inasmuch as the coordinate axes rotate during 
the motion, the time derivatives of these quantities do not represent 
the actual airplane accelerations, which are required in analyses of the 
loads experienced by the airplane and the degree of passenger discomfort. 
If the deviations from a mean flight path are assumed to be small, the 
accelerations normal and parallel to the chord or longitudinal axis of 
the airplane are substantially the same as the absolute vertical and 
horizontal accelerations z and x, which can be obtained from the 
relations 

z = -Wp + U§ 
x = -Up - g0 

The transfer functions for these quantities can then be obtained by 
introducing these relations into equation (43). 

In studies of the longitudinal stability of airplanes, equation (43) 
Is rarely solved in the form given here. It is usually reduced to two 
equations with two unknowns, either u p and 0 (the phugoid case) or 

Wp and 9 (the short-period case), the short-period case being usually 

the one of primary interest. The part of the turbulent energy contained 
in the frequency range near the phugoid frequency is relatively small, 
so that the phugoid case has no significance for the analysis of loads 
and accelerations resulting from atmospheric turbulence. Hence, the 
short-period case, which ignores the phugoid oscillations, furnishes an 
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excellent approximation to the loads and accelerations associated with 
the longitudinal motions of an airplane in turbulent air. However, 
another two-degree-of -freedom, case, the one involving Vp and Up, is 
useful in certain studies of the effects related to the interaction of 
horizontal and vertical components of turbulence. 

Both of these two-degree cases can be reduced to the s ingle -degree- 
of -freedom case involving only z (or Wp) . For airplanes which have 

a laxge moment of inertia in pitch this simple case furnishes a good 
approx ima tion - It has been studied in reference 2, where substantially 
the same approximations to the unsteady-lift effects were made as are 
made here, except that in reference 2 apparent-mass effects were included 
(These effects are not included in the stability derivatives used in 
equation ( 4 - 3 ) because they are usually small - less than 1 percent of 
the mass of the airplane - and are different for each degree of freedom. 
However, if desired, the apparent mass pertaining to a given degree of 
freedom can easily be added to the airplane mass in calculating the 
stability derivatives . ) However, this approximation iB more nearly valid 
for calculating peak loads persuant to an entry into a sharp-edge gust 
than for calculating the response due to random turbulence. Consequently, 
in the following sections, attention will be confined to the short-period 
case, although the analysis is equally applicable to the other case and 
easily extended to the case of three degrees of freedom. 


Solution of the Equations of Longitudinal Motion 


Transfer functions .- If the degree of freedom pertaining to x 
(or Up) is ignored, the solution of equation ( 43 ) can be written as 


z 

► = 

Bgfro) 

Ha(®) 

■ 

* * 
w 

.1 


lf(®) 

B^(o>) 


u 


where the transfer function Hg(cu) is defined by 


= fl(k) 


, w p A „v 
-Agm^ + B”is> + C‘g 

-A + BqJo + Cq 


( 44 ) 


( 45 ) 


where, in turn, the coefficients are defined in terms of the stability 
derivatives (see table l) by 
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1 — 1 
II 

«? 

> 

I 

II 

B 0 = "(Zw + M q + I3M*) 

Bf = Zy'(U% + 1%) 

Cq - MgZy - UMy 

Qg = -U(My>Zy - MyZy 1 ) 

V/ \ 

The transfer function can be defined similarly in terms 

coefficients 

> 

it 

+ Zy'Mtf 


% = -(V^w “ ^w^w') 
"vr 

% - 0 


In these equations a distinction has "been made between the values of Zy 
and My- which occur on the right side of equation (43) and are here 
designated by a prime mark, and those on the left side of that equation. 
The primed derivatives pertain to the lift and moment directly due to 
gusts, whereas the unprimed derivatives pertain to the lift and moment 
due to airplane motion; the reason for this distinction is discussed in 

u u u u u 

a later section. Furthermore, the coefficients A”, Bg, Cy, Ag, Bg , 
U. "W "VT 

and Cg“ are the same as the coefficients A”, B", . . . except that 

Zy, Zy’, My, and My' are replaced by Z^ and M u r . 

With these transfer functions the mean-square values of z and 
can be calculated from the spectra of w and u by using equation (l), 
provided that the simultaneous action of w and u is taken into 
account. In order to analyze this effect the vertical acceleration z 
will be considered, but the analysis will be applicable to 9 or any 
other characteristic of the airplane which responds to w and u. 
Furthermore, the transfer functions need not be those considered in 
equation (44), but can be those calculated for the three-degree-of- 
freedom system or for a flexible airplane. 
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Fourier transforms of the transfer functions (see eq. (3a) ), are more 
convenient. In terms of these Indicial-response functions, 

p CO p 00 

z(t) = / ti|(t 1 )w(u(t-t 1 ) ) dt L + / h" (t-L) u^Uft-tj) dt L 

U —CO W —00 

Then, if w(t) and u(t) are stationary in a statistical sense, the 
correlation function for 2(t) can he written as 

V T > =/ f SOt) ^(*2) +v( U ( T+t l“ t 2)) + ^(tl) ^ 2 ) t u (U( T -rtl-t2)) + 

—CO —CO 

Hg(ti) bi(t 2 ) ^u(TI( T +t L -t 2 )) + h»( tl ) hl(t 2 ) tvu(u(-T+ti-t 2 ))] at! dtg (46) 
where ^^(Ut ) is a cross correlation of w and u defined hy 


t w (Ur) = w(x-KJt) u(x+u( t+T ) ) 

Now, if the turbulence Is isotropic, the mutually perpendicular 
velocity components u and w at points in the XY -plane are statis- 
tically Independent, and their cross correlation is zero. Therefore, 
the two terms in equation (4-6) Involving ^^(Ut) vanish, and the power 

spectrum of z Is 


<Pg(<o) = |^(“)| 2| ? v (“) + |a|(a))| 2 cp u (ai) 


(47) 


so that, generally speaking, the power spectrum of a response which 
depends on both the horizontal and the vertical component of turbulence 
Is simply the sum of the power spectra of the two contributions, provided 
the turbulence Is isotropic. (This statement can be shown to be true 
even If the distribution of the gusts over the span is taken into account.) 


For the short-period two-degree-of -freedom case, then, the contribution 
due to w is 



i*«i 2 


M3 

! CD — 

- ( 

$ 2 ] 

» 2 * (05 f 

ID* - 1 

V 2c 0 - B 0 2 ’ 


1- c 0 2 


<P„.(cd) do 


and the contribution due to u can be obtained from the same expression, 
but with the subscript and superscript w replaced by u. However, the 
ratio of the two contributions Is in the order of 4a, 2 , where a Is the 
trim angle of attack in radians, measured from the zero-lift conditions. 
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Consequently, except at very high lift coefficients, such as those used 
in landing, the contribution due to u ; is usually negligible compared 
with the one due to w and is disregarded in subsequent sections. 

By means of the lift-attentuation function given in equation ( 7 ) 
and the point spectrum used in the first part of the paper, the Inte- 
gral for z 2 has been evaluated (by using the technique of partial 
fractions for the integrand) for the example of reference 11. (The 
lift -attenuation function of equation ( 7 ) has been used despite its 
shortcomings In order to facilitate the analytic integration of the 
spectra. ) The results are shown in figure 6, as are the results calcu- 
lated similarly for the three-degree-of -freedom case, the other two- 
degree-of -freedom case (horizontal and vertical motion; referred to In 
the figure as the zero-pitch approximation), and the single-degree-of- 
freedom case (vertical motion) . An examination of this figure Indicates 
that, at least for this airplane, inclusion of horizontal motions does 
not affect the mean-square normal-acceleration response to any significant 
extent... These calculations pertain to a "small" airplane, and no dis- 
tinction has been made between the primed and unprimed values of Zy 
and My. 


The preceding treatment of the short-period caBe haB the advantage 
of using readily available information concerning the characteristics of 
any given airplane. For the purpose of trend studies a dimensionless 
form of the transfer functions is preferable. The preceding equation 

for z 2 can be written in dimensionless form (the contribution of 
horizontal gusts again being neglected) as 




w 2 _4_ 

U 2 k 2 


riifoof 2 

Jo 


+ 4( 

V -1' 
IW 


CVJ 

1 

-Si 

(ko 2 - V 2 

)k? + l 

IwW) 2 


u cp v (co) 


dk 1 


(48a) 


and, similarly. 


(2gA0 2 






_j£ u J^ldk» 

v 2 )k 2 + (k Q 2 + v 2 ) L#w2 

(48b) 


where k Is the mass parameter 


_ 8m 



N 


NACA TN 3910 


47 


v is the d im ensionless damping coefficient 

v - 0 lQ Se 2 
-2U T l/2 

is the time to damp to one-half amplitude, which is given by 
log e 2 


so that 


= ■ IK + ^ + ™s) 

v.lfx 1 me 2 % + 

“ V ' 2 ly c in, / 


and ko is the dimensionless frequency of the short-period oscillations, 

one 

110 = W~ 

with 


oIq - ^MqZ^ - UMff - + Mq + UMw)* 


so that 



,^ l ^ +2 s 

x y \ c iu c i<t. 


1 mc^ + C m d 

2 % C rc 


2 


Thus, for this two-degree-of -freedom case the dimensionless mean- 
square responses are functions of only two additional parameters, which 
are dimensionless forms of the main characteristics of the short-period 
case (the short-period frequency and the time to damp to one-half ampli- 
tude), beyond those encountered in the single-degree-of -freedom case, 

r* 

namely, the mass parameter k. and the scale parameter s = 


Bending -moment spectrum and mean-square value .- For a small airplane 
the instantaneous bending moments at various points on the span are 
proportional to the instantaneous normal acceleration. For instance, 
the root bending moment can be written as 




(k-9&) 
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where and % are, respectively, the lift on and the mass of the 

wing, and where y and y are the lateral distances to the center of 
pressure of the lift on one wing and the center of gravity of the mass 
of one wing. However, is proportional to the lift L on the entire 

airplane, which in turn is proportional to^the normal acceleration, so 
that 

Ly. = mz 

and 

% = §(^ ? - j y )] 2 

Hence, the spectrum and mean-square value of Mg are proportional to the 
spectrum and mean-square value of "z, respectively, the constant of 
proportionality being the square of the quantity in the brackets of 
equation (49b ) . 

Special Problems Related to the Longitudinal Motion 
of Large Airplanes 

Single-degree -of -freedom, case .- In the preceding sections the air- 
plane has been assumed to be small in the sense of this report, and 
neither the Instantaneous lateral variation of the intensity of turbulence 
nor the difference between the Instantaneous intensities at the wing and 
the tail has been taken into account. In this section this restriction 
is removed by introducing the aerodynamic forces calculated in the first 
part of this report into the dynamic analysis of the preceding sections. 
The arguments advanced in the preceding, section for ignoring horizontal 
gusts and horizontal motions are equally valid for the large airplane; 
therefore these gusts and motions will not be considered here. 

For the single-degree-of -freedom case involving only vertical motion 
the required modification for the normal acceleration is very simple. 

For this case the transfer function is 

I&o) = -fS(k) JlL— (50) 

16 icu - Zy 

where attenuated by 0(k ), is the stability derivative for vertical 

gusts and, hence, represents the lift per unit gust intensity. Therefore, 
if the result for the lift calculated In the first part of this report 
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is used, the mean-square normal acceleration "becomes 


0 


ol 


2 .2 


2 Zv®' 


05 2 + Zy - 2 


cp„. (oi) do 


(The function cp We (to) is defined in the first part of this report.) 

This expression differs from, the result obtained in reference 2 only in 
that cp-yCoO is here replaced by cpy e (aj) . 


However, even for the single-degree case the calculation of the 
bending moment now becomes a considerably more complicated problem, 
because the lateral centers of pressure of the lifts due to the motion 
of the airplane and directly due to turbulence no longer coincide. 
Equation (49a) now becomes 


Mg = 


— T Ct Q.S • 

£ fw V z 
2 L U U 



+ 



(51) 


where is the lift-curve slope for the entire airplane, h L (t) is 

the response function used in equation ( 38 ) , and / M (y) is the function 

M(y,q) used in equation ( 38 ), with y = 0. In the first two terms on 
the right side of equation (5l), z and z can be related to w by 

TnpflnR of h” (the Fourier transform of the function Hg presented 
in eq. ( 50 )) and the lift directly due to turbulence obtained in the first 
part of this paper, so that these two terms can be written in the form 

r°° r>b /2 

bj_«, ^ J b / 2 7(y) v ( U ( t " t 0' y ) <3y 

and, hence, equation ( 51 ) becomes 


M B (t) = dt L J ^ jhL(t L ) 7 M (y) + ^(* 1 ) 7(y)] w(u(t-t L ),y) dy 
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The required influence function for the bending moment is, thus, 


^(y^t) = J[h L (t) 7 M (y) + h^t) /(y)] 


( 55 ) 


Although the two terms of this function: are products of a time-dependent 
and a space -dependent constituent, as in equation (5), their sum cannot 
be split up in this manner . Hence, the techniques used in the first 
part of this paper are not directly applicable. The general treatment 
of problems involving influence functions for which the assumption of 
equation (5) is not valid will be considered in the last part of this 
paper. However, in this section a special technique will be used that 
applies to cases for which the influence function can be expressed as a 
sum of several terms (two in this case), each of which can be expressed 
as a product of two functions, which depend, respectively, on time alone 
and distance along the span alone. (See also ref. 13 •) This approach 
is more convenient than the general approach of the last part of this 
report when the number of terms is two or, possibly, three. (Although 
in this section only two terms will be considered, the generalization 
to three or more terms is straightforward.) When the number of terms 
is greater than three, the general approach becomes more convenient. 

If the value of M B given by equation (52) for t + t iB multi- 
plied by the value for t and the result is averaged, and if the Fourier 
transform of the resultant correlation function is then taken, the fol- 
lowing expression is obtained for the spectrum of the bending moment: 

c Pm b ( cu ) = |H L (^)| 2 9 W e 2 (a)) + !%(“) |%r e ^ + 211 %*(“>)} (5^) 


where the symbols R 


{} 


and 


designate, respectively, the real part 


and the complex conjugate of a complex number. In this equation %(cu) 
and H b (cd) are the Fourier transforms of -hjyj(t) and h L (t), c P V r e (cj) 
is the previously defined averaged spectrum for the lift, and cp^ 


and qv e2 
the value of 


are averaged spectra obtained in a similar manner but with 
T given by equation (12) replaced by 


^(ri) 

r 2 U) 



L 


(b/2)-T] 


L 


■b/2 

(b/2)-rj 

■b/2 


7 ( 7 ) 7 M B (y +T l) 

[7M s (y) 7M s (y +T i) + 7Mjy> ^(y-^)] 


■\ 



(55) 


where the subscripts s and a refer to the symmetric and antisymmetric 
parts of 7 m , respectively. 
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The contribution of the antisymmetric part of 7 ^ to r 2 and, 

hence, to the spectrum for the bending moment stems basically from the 
asymmetry of the instantaneous distribution of gust intensity over the 
span. This asymmetry gives rise to a rolling moment (which was consid- 
ered in the first part of this report) and, hence, to rolling motions, 
which contribute additional bending moments due to the aerodynamic and 
inertia loads associated with these motions. If, for the purpose of 
calculating the bending moment due to symmetric flight through turbulent 
air alone, these motions are disregarded, then the contribution of 7 ^ 

to r 2 should be disregarded as well. The problem of combined symmetric 
and antisymmetric motion will be considered presently. 

As pointed out in reference 13, the mean-square bending moment calcu- 
lated in this manner may be smaller or greater than the value calculated 
by ignoring spanwise variations in gust intensity. If the mass of the 
airplane is almost entirely contained in the fuselage, the decrease in 
the lift which results from taking these variations into account causes 
a decrease in bending moment. However, if most of the mass is in the 
wing, the net bending moments (aerodynamic less inertia) for a uniform 
spanwise gust are very small, and the effect of taking spanwise variations 
of gust intensity into account is to increase the mean-square bending 
moments . 

Two-degree-of -freedom case .- For the two -degree-of -freedom (short- 
period) case the analysis given for the normal acceleration in the pre- 
ceding section can be extended as follows: As indicated in equation (41) 

the transfer function for z is now 

Bg(o>) = 


(ihe following analysis can be applied equally well to 9 by using 
E^(cjo) instead of H^(o)).j 


! cd 2 + (UMfr + Mqjirn + UmJ 


m p 

-a i + B 0 ia) + C 0 


m(-V) 0(k) 


U(-Zw) 


% -CD 2 + BqIco + Cq 


IyV 0 (k) 


( 56 ) 


If the expressions inside the brackets of equation ( 56 ) are desig- 
nated, respectively, by Hp(cu) and H 2 (cd), and their Fourier transforms 
by h-j_(t) and h 2 (t), then, as a result of the definitions of Z^’ 
and My.’ , 


/ , 00 

hpOa) L(t-t-,} dt x 

-CO / 



hg^i) M(t-t x ) dt-L 


( 57 ) 
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where L(t) and M(t) are the instantaneous lift and pitching moment 

due to the vertical component of atmospheric turbulence , which have been 

obtained in the first part of this report. The calculation of 9 .. (as) 

z 

Tig 

or of z thus requires not only the spectra of L and M, the calcu- 
lation of which has been discussed, but also the cross spectrum of L 
and M, which has to be calculated directly from equations (4) and (42a). 
The result is 


= |H L (a3) j 2 |h(cd) j 2 9„. e (aj) + | 2 jjB^Cco) J 2 tJV e (aj) + 

| j 2 cp w .(ao) + 2 R je ^ 13 E^cd) H w .*(n)| 9w e # (^0) 
2|r|h 1 (o>) E^oj) H(a>) H^*(a>) j qv e (a>) + 

RjHxCcu) Er>*(ca) H(ai) j qv e *(^0)J 


(58) 


where the first two terms represent the contributions of the spectra 
of L and M, respectively (see eqs. (l4) and (42b)), and the third 
represents the contribution of the cross- spectrum of L and M. For 
the present purpose the functions E(cs), E^cd), and E^cd) can be 
expressed as 

H(o) = m(-Z w -) 0(k) 

Hy-tn) = m Ax a (-£„•) 0(k) 

Ht(o>) = [lyyJV " m ^a (-Zv)] 0(h) 

where Ax a is the distance from the aerodynamic center (tail off) to 

the airplane center of gravity. This definition of the contributions of 
the wing and tail to the pitching moment is based on the considerations 
that the direct contribution of the wing can be estimated with good 
accuracy and the total pitching moment is likely to be known from experi- 
ments, so that the contribution of the tail (which includes the effect 
of the wing lift on the downwash at the tail) can be determined as the 
difference of the two. The functions 9 ^( 0 ) and 9 w - e *(o,0) have been 

defined in the first part of this report. In view of the fact that the 
function 0(k) contained in same of the terms of equation ( 58 ) always 
appears in terms multiplied by others which contain 0*(k), only the 
absolute square of this function is required, as before. 
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Special Problems Related to the Lateral Motion of 
Large Airplanes 

Hie equations of motion in the lateral degrees of freedom (roll, 
yaw, sideslip) have the same form and can be solved in the same way as 
the equations for the longitudinal motion. (See pp. 111-53 to HI -67 
of ref. 11.) Again it is convenient to cast the problem in the form 
used in a stability analysis in order to take advantage of the results 
of such an analysis. For a small airplane it is necessary only to 
replace the terms due to rudder deflection by corresponding terms involving 
side gusts, namely, 

lg R S R by Y v §0’(k) 

^r s R by 

% r 8r Ly N p J0’(k) 

where 0 1 (k) is a suitable side-force attenuation function for side 
gusts and the notation of reference 11 is used for the other terms. 

The terms corresponding to aileron deflection can be disregarded for 
the small airplane. In the lateral degrees of freedom the small airplane 
thus reacts only to side gusts. On the other hand, the large airplane 
also reacts in the lateral degrees of freedom to vertical gusts through 
the rolling moment calculated in the first part of the present paper. 

If this rolling moment is to be included, it replaces the term 6^ 

used in reference 11. 

Instead of treating all three degrees of freedom simultaneously, 
in stability analyses two one-degree-of -freedom cases are often consid- 
ered, namely, the one of sideslip alone, with angle of yaw equal and 
opposite to angle of sideslip (the Dutch roll case), and the one of 
rolling alone. 

The Dutch roll case may be used for gust-load purposes in connection 
with yawing and sideslipping motion due to lateral gusts, provided the 
phase of the motion is not important. (As may be noted from the preceding 
sections, the phase of a transfer function is important only in terms 
Involving cross spectra. ) Also, the Dutch roll case may prove useful 
in calculating the vertical-tail loads resulting from flight through 
turbulent air, particularly if the human or automatic pilot holds the 
wings substantially level by means of the ailerons. However, In general 
it does not appear to be as satisfactory an approximation as the one- 
and two-degree-of -freedom approximations for longitudinal motion. 
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For rolling motions due to rolling moments , the single-degree case 
of rolling alone appears to furnish a very good approximation. Although 
the rolling motion causes yawing and sideslipping motions, these motions 
do not appear to reflect on the rolling motion. Thus, the rolling motion 
which results from the rolling moment can probably be calculated fairly 
accurately without regard to the other lateral degrees of freedom. 
Furthermore, within the assumption of small motions, the stresses asso- 
ciated with these other lateral degrees of .freedom do not generally 
contribute appreciably to those associated with the longitudinal degrees 
of freedom in the parts of the structure for which the latter are critical, 
such as the wing (although they may be critical for other parts of the 
structure, such as the vertical tail). Therefore these degrees of freedom 
(yaw and sideslip) will be ignored in the treatment of the large flexible 
airplane in the last part of this paper. However, if chordwise bending 
effects (deformations parallel to the chord) are important, as they may 
be in some cases at speeds close to the flutter speed, these other lateral 
degrees of freedom may have to be included in the analysis. 

For the large airplane, which responds in the lateral degrees of 
freedom both to vertical gusts and to Bide gusts, the superposition of 
the resulting responses, such as stresses, may be effected in the way 
indicated for the interaction of horizontal and vertical gusts. If the 
turbulence is isotropic, the vertical and lateral gusts are statistically 
independent for points in the XY-plane, so that the spectrum of a given 
response is equal to the sum of the spectrum of that part of the given 
response which is due to vertical gustg , and the spectrum of that part of 
the response which is due to lateral gusts , 


Combination of the Results Obtained From Analyses of the 
Longitudinal and Lateral Degrees of Freedom 

The instantaneous wing stresses depend both on the motions in the 
longitudinal degrees of freedom (primarily vertical motion and pitching) 
and on those in the lateral degrees of freedom (primarily rolling) . The 
purpose of this section is to indicate how the stresses associated with 
vertical motion and pitching can be combined with those associated with 
rolling due to vertical gusts, particularly in the case of a large rigid 
airplane. (A small airplane, flexible or rigid, does not roll as a 
result of the action of vertical gusts, and for the large flexible air- 
plane it is more convenient to consider rolling motion simultaneously 
with the other motions, so that the superposition is effected automatically 
in the process of ..-obtaining the required transfer functions.) For all 
airplanes the effect of side gusts can then be taken into account, if 
isotropy is assumed, by adding the stress spectra directly. 
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In this section the instantaneous stress at a given part of the 
wing will he assumed to he proportional to the instantaneous bending 
moment at that section of the wing, so that consideration can he confined 
to this bending mamentj also, the airplane will he considered to he free 
to move in only two degrees of freedom, namely, vertical motion and 
rolling. The extension of the following arguments to stresses which 
depend on the vertical shear and the torque as well, and the inclusion 
of pitching as an additional degree of freedom, can he effected readily 
and will not alter the conclusion reached here. 

When rolling motions of the airplane are taken into account, the 
bending moment due to the aerodynamic forces associated with rate of roll 
and to the inertia load associated with rolling acceleration can he 
expressed in terms of the rolling moment by a superposition integral. 

In turn, the rolling moment can he expressed in terms of a superposition 
integral involving the instantaneous gust intensities along the span, 
the influence function y'(y) used in the first part of the paper, and 
an associated response function b^(t'). Hence, equation ( 52 ) is modi- 
fied by the inclusion of a third term and becomes 


*b/2 


MgCt) mj at-L J £h L (ti) 7 M (y) + 7 (y) + l^(t L ) 7 , (y)]v(u(t-t 1 ),y) <3y 


If the spectrum for Mg is now calculated in the manner used pre- 
viously, the following result is obtained: 

9^(01) = |H L (o))| 2 q> We ^(o 3 ) + |%(o>)| 2 9 We (a>) + 2R |h l (oj) %*(a>) J (a>) + 

|^j(“)| 2< Pw e3 (“) + 2R |h l (o)) q^(<o) (59) 

where B^(cs) is the Fourier transform of hjJj(t), and where cp w (co) 


and (o>) can he obtained from equation ( 12 ) with the following values 

of and r^, respectively: 




(h /2 ) — T] 


Vi> " lf_ 


■h/2 
(b/2)-7] 
■h /2 


7*(y) 7’(y+Ti) <3y 


7’ (y) 7 Ma (y +T i) ay 
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(The function I\j(t)) is four times the function T r ( T| ) considered 

previously in connection with the rolling moment. This factor of four 
must be taken into account in l^(t).) 


The spectrum cp^ may be considered to consist of two parts: The 

terms listed in the first line of equation ( 59 ) , but excluding the con- 
tribution of 7 ^ to (p w __ (see eq.. ( 55 )), represent the contribution 


T v ( 




of the symmetric parts of the instantaneous gust distributions or the 
contribution associated with symmetric motion; the terms on the second 
line and the contribution of y u to cp^ represent the contribution 


'“a 


e 2 


of the antisymmetric parts of the instantaneous gust distributions or 
the contribution associated with rolling motion. Therefore, the power 
spectrum of the stress due to gusts and combined symmetric and rolling 
motion resulting from the gusts is the sum of the two power spectra 
(that for the gusts and symmetric motion alone and that for the gusts 
and rolling motion alone), provided the direct contribution of the gusts 
is split up into a symmetric and antisymmetric part and each is taken 
into account only once, in connection with the appropriate type of motion. 
The cross -correlation terms between the symmetric and antisymmetric con- 
tributions to the stress can be shown to involve integrals which contain 
products of symmetric and antisymmetric influence functions and, hence, 
vanish, so that the cross correlations are zero; hence, the two parts 
of the combined spectrum are statistically independent and, therefore, 
directly additive. 


DYNAMICS OF THE SMALL FLEXIBLE AIKPLANE 


The purpose of this part of the paper is to consider the transfer 
functions relating the stresses at various points of a small flexible 
airplane to the vertical gusts which cause them. The longitudinal and 
lateral degrees of freedom are still separable, and only the longitudinal 
degrees will be considered; the lateral degrees, which are involved in 
analysis of the response to side gusts acting on the vertical tail, can 
be analyzed in the same way. Therefore, for this case only one-half of 
the wing need be . considered as a result of the symmetry (or antisymmetry, 
in the case of the lateral degrees of freedom) of the problem. 

The method which is outlined in this part consists in an extension 
of the numerical-integration method of aeroelastic analysis described 
in reference lit- to sinusoidal motions of the airplane. This extension 
takes into account the facts that the aerodynamic forces now have out- 
of -phase. as well as in-phase parts and that vertical deflections must 
now be calculated separately because the structural deformations can 
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no longer be characterized, by angle -of -attack changes. Also, the rigid- 
body" degrees of freedom (vertical and pitching motion of the airplane 
as a whole and structural deformations of the tail) are now taken into 
account. 

The result is a set. of linear algebraic equations (which serve as 
an appro xima tion to the ordinary differential equations that characterize 
the problem) for the airplane motions and deformations in terms of the 
applied aerodynamic forces directly due to guBts. The desired transfer 
functions can then be obtained from solutions of these equations at 
various frequencies , and the power spectra of the stresses are given by 
the product of the absolute square of these transfer functions and the 
point power spectrum of the vertical component of turbulence. The same 
equations may be used to calculate with little additional effort certain 
aeroelastic effects, such as the static aeroelastic deformations and the 
flutter speed, which are usually obtained in separate analyses. 


Loads Applied to the Wing 

The loads applied to the wing stem from three sources: The aero- 

dynamic loads directly due to the action of the gusts, the aerodynamic 
loads due to the motions of the airplane, and the inertia loads. 

The lift and pitching moment (about the elastic axis) per unit span 
on a two-dimensional airfoil undergoing sinusoidal angle- of -attack changes 
and vertical motions in incompressible flow are (see ref. 15 ) 


Za = acqc jc(k) [(1 + 2^ik)cc - Ik + (*£ - e^a + f ^ j 

“a = 2 *qP 2 jcOO e ^ 1 + 2 % Ik)a - ** + ^Y^ “ J ^ 


(60) 


The terms multiplied by C(k) are referred to as the circulatory terms 
because they are calculated from the bound and shed vorticity, and the 
others are referred to as the potential terms . The potential terms are 
in the na ture of additional-apparent-mass effects , and all those that 
involve k 2 are usually treated together with the inertia forces rather 
than with the aerodynamic forces. For compressible flow, however, the 
forces are calculated in a different manner, and the division of the 
forces in to circulatory and potential parts then has little meaning. 
Consequently, in order to facilitate the extension of this analysis to 
compressible flow, thiB distinction will not be made herein. 
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The aerodynamic forces are therefore -written as 


Z a = 2itq.cj& 1 (k) a + C 2 (k) -yy| 
1 % = 2jrqc 2 j8' 5 (k) a + C^(k) 


so that for incompressible flow 

C x (k) = (l + 2e^ik) C(k) + - e^k 2 

? 2 (k) = -ik C(k) + ^ 


^(k) = (1 + 2e^ik) e x C(k) - ^ ik + + e^jk 2 

C^(k) = -ike-j^ C(k) - k 2 


(61) 


In order to calculate the lift at a given point of a wing of finite 
span an appropriate Green's function is required. An approximation to 
this function based on a reciprocity theorem of linearized lifting-surface 
theory. Is given In the first part of this paper j as used for the compu- 
tations of that part, this function implies the assumptions that the 
spanwise distribution of the lift for oscillations of the wing as a whole 
Is substantially Invariant with frequency and that this distribution is 
the same in direct as in reverse flow. Neither of these assumptions is 
essential to the analysis but both, and particularly the first, simplify 
It considerably. With these assumptions, the desired lift distribution 
is then given by expressions of the form of equations (32a) and (34). 

For the present purpose, however, a set of aerodynamic influence 
coefficients is required, rather than influence functions. Such a set 
of coefficients, based on the same ideas, can be obtained readily by the 
techniques used in references 9 and 10. The result may be expressed as 
follows: - - 

f . 

{ l }a = C L (k)[q]{a} + C 2 (k)[q]|^-| > ( 62 a) 

where the aerodynamic-influence-coefficient matrix [q] is defined by 

M = (1 - K)t7l{l}LUl/l + K 
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■where, in turn, is a unit column matrix, and |_I J is a row of 

integrating coefficients suitable for integrating a continuous function 
for a range of its argument from 0 to 1. Urns, for instance, if n 
equidistant points on the semispan are considered, and n is odd, then 
according to Simpson’s rule 


LiJ 


1 1 k 2 k 

n - l[3' 3' 3 J 3' ' ‘ 


t if 

y 3j 


Very little information is available concerning the spanwise distri- 
butions of the pitching moment on wings of finite span in unsteady flow. 
By means of the reciprocity theorems an appropriate Green’s function 
could be estimated if the lift distribution for wings with parabolic . 
camber were known; however, such lift distributions do not appear to have 
been calculated for wings of finite span in unsteady flow. In fact, 
relatively little is known about moment distributions even in steady flow. 
However, the available information indicates that the local center of 
pressure does not appear to be very sensitive to the lift distribution. 
(See ref. 9, for instance.) It will be assumed that this is also true 
in unsteady flow at a given frequency, and that, furthermore, these cen- 
ters of pressure are given by two-dimensional theory. With this assump- 
tion the moment can be written as 

{™T} a = Cl^ccL jc 3 (k)tc} [Q]{a} + fy(k) ]>} [ q] j ( 62 b ) 


The lift and moment distribution due to the gust can be calculated 
in a similar manner. In the following analysis, the magnitude of the 
gust intensity Is considered to be unity, and the longitudinal reference 
point is the intersection of the elastic axis and the wing root, so that 
the instantaneous gust intensity at any station y is 

-ikM tan A 

w = e c '^ ( 63 ) 


With this function w. 




Cl^cq 0(k) 


U 


MW 




(64) 


where 0(k) is the Sears function, as before. 
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Finally, the inertia loads are 

2^ = -Hz + m^ca 


m T i = 



( e 2°) 2 ) 


mcc 


/V 

me. 


,cz 



The loads applied to the wing can thus he written, in summary, as 


i 


20^ %>(*) W + 
2\2 / 

c L (k) [q] - 
- k - 2 - Tme 2 c~| 

m 

I*l 4 

80 

W hi [Q] - 

Cr^c C 5 (k) hi [Q] -l- 

f 2 I r2 + (vflR 
id) 

w 


[megc] 

m 


u L L _r J 


■SW i 


( 66 ) 


Loads Applied to the Tail 

The loads applied to the tail are similar in nature to those applied 
to the wing, hut the tail experiences additional loads as a result of 
the downwash produced hy the lift on the wing. Again, little is known 
about the downwash in unsteady flow, and even in steady flow the downwash 
cannot he predicted accurately because of houndary-layer effects on the 
fuselage and the wing root. Consequently, even in steady-flcw analyses 
experimental results are usually relied upon. 

In the following analysis, the assumption is therefore made that 
experimental results are available for steady flow, in the form of the 
downwash derivative de/cfcc,. In order to determine the attenuation of 
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this value with frequency, the results of the analysis of reference l£ 
will he used. These results indicate that the time variation of the 
tail lift due to the downwash caused hy the wing lift which results from 
a unit .Jump in the wing angle of attack can he approximated hy an imme- 
diate Jump in the tail lift of -0.l6 of the steady-state value and another 
jump to the steady-state value after the time required to travel the dis- 
tance from the 45 -percent-chord point of the wing to the quarter-chord 
point of the tail plus another eighth of the chord length. Hence, for 
sinusoidal angle-of -attack changes the tail lift due to downwash is 



a,z 


-C 





(l+ ae^ikja,,- ik 



(6?) 


where x-^' is the distance from the intersection of the elastic axis 
and the wing root (assumed for this purpose to he at the percent point 
on the root chord) to the aerodynamic center of the tail. As pointed 
out in reference 16 , this approximation is valid only for k < 0 - 35 ; 
however, this range is adequate for the present purpose. 


Similarly, the downwash associated with the wing lift due to sinus- 
oidal gusts gives rise to a tail lift which, within this approximation, 
is 


K) g - £ | 


-ik 


- 0 .l 6 + l.l£e 



- 0 . 6 ik 

c 


w„ 


( 68 ) 


(The additional lag represents the time required to travel the 0.6 root 
semichords from the 45-pereent-chord point of the root, which is the 
reference point for the gusts, to the 75-percent-chord point of the wing 
root, which is assumed to he the point governing the lift at the w ing 
root, inasmuch as it is the centroid of the influence function for the 
chordwise pressure distribution. ) 

The other aerodynamic forces are those due to the motions of the 
airplane, those due to the tail deformations, and those directly due to 
the gusts. On the basis of the assumptions made in the preceding section, 
these forces are 

^a = G Ip^t s t[%( k > + %(h) + <* 7 ( k ) Ax + c 6 (k) ( 69 ) 
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where 

( 1 c t\ t c t 

1 + 2ik — §-) cfe k) + t _ k 2 

^ c / \c / 2 c c c 

C 6 (k) = -ik k) + J ^ k 2 

S T ( k) . (x + ** 3£) 0^ *) * | S * - e 3t (i) 2 ** 

and Oj, and Zj, are the angle of attack and vertical displacement of 
the airplane at the wing root. -- 

The inertia load on the tail is 

^i = “ m t (k"r + Az - xt’c^) 


or 



(71) 


Here the center of gravity of the tail has been assumed to coincide with 
its aerodynamic center; in order to remove this assumption it is necessary 
only to add (or subtract) the distance between the two to x^.' in equa- 
tion (71) • 


The normal forces on the tail can then be summarized as follows: 
r*t = q.[ji(k) z r + F 2 (k) + ^(k) Ax + F^(k) Az + F^(k) w t J (72) 


where 


<lt- s t 
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F 2 (k) ‘ “lot T S H ef 5 (k) ' ( k + £ v k ) jg 
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The pitching moments corresponding to these normal forces can he 
obtained in the manner employed for the wing. However, inasmuch as the 
tail chord is usually small compared with the fuselage length, the 
travel of the center of pressure of the tail is small compared with the 
length x^'. Hence, the center of pressure will be assumed to remain 
at the aerodynamic center of the tail, and the pitching moments are 
then -x^' times the corresponding forces, so that 

M t = “^xt' [ F i( k ) z r + F 2 ^ k ) “r + F 3 ( k ) ^ + F lj.( k ) + F 5 ( k ) (73) 


Wing and Tail Deformations 

The wing deformations may be calculated either from structural 
influence coefficients or from the bending and torsion stiffnesses of 
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the wing used in conjunction with simple beam theory. The latter approach 
will be followed here, based on the method of reference 12. 


The bending and torsion moments on the wing structure may be obtained 
by integrating the applied loads. If numerical methods are employed to 
perform these integrations the results may be written as follows: 



where the integrating matrices (see ref. 12) perform the following 
operations : 

[i]{f} o J ^ f(i') as' " 

Lll] {f } •£= / 7 " f(l") di" di' 

[ 1 "] jfj =0 J* f(g') dg’ 

[n"]-jf| o f (s B ) dg" ag' 
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These integrating matrices may be based on the trapezoidal rule, Simpson's 
rule, or any other numerical method; the intervals chosen for | 

(0 ^ |=l) need not be of constant width unless a specific rule demands 
a uniform spacing. 

The structural deformations of the wing may then be written in 
terms of the applied loads as 



where the submatrices j(T)] , [©]•> [d)] > [®] > [@] > and [©] 

designate, respectively, the four quadrants and two halves of the matrices, 
obtained by postmultip lying the square matrix of equation (75) by the 
square and the rectangular matrix of equation (66), respectively. 

For the purpose of the present analysis, which is concerned primarily 
with the wing stresses, the tail deformations are treated by including 
only the vertical displacement and angle -of -at tack change of the tail as 
a whole due to the tail load. These quantities may be obtained from a 
static test which consists in applying a concentrated normal load at the 
aerodynamic center of the tail and measuring these deformations. They 
may also be obtained from a vibration test in which the deflection curve 
of the rear part of the fuselage in the lowest vertical-bending mode is 
measured; in this case the desired spring constants can be deduced from 
the frequency relation of a simple mass oscillator in terms of the meas- 
ured frequency and of the mass of the empennage (including the part of 
the fuselage which may be considered to move with the empennage) . In 
the absence of such tests these constants may be calculated in an anal- 
ogous manner. 

The tail deformations may then be written as 

Ax = -KjLt 
Az = 




( 78 ) 


so that, also. 



66 


MCA TN 5910 


Inasmuch as these deformations are not independent of each other, only 
one need he retained in the analysis. Therefore, if Az is eliminated 
by means of equation (78) , Ax can then be obtained in a form similar 
to that vised for the wing deformations in equation (78): 


Ax = -qKq F-]_(k) z r + F 2 (k) °S: + 


F,(k) - 


Ko 

£ 


Ax^ - qK^ F s (k) w t 
(79) 


Equations of Motion 


Equations (76) and (79) are equations of motion inasmuch as they 
describe balances of aerodynamic, structural, and inertia forces. In 
fact, if the airplane fuselage were immobile (a^ = z r = 0) , they would 

be sufficient to calculate all unknown quantities. However, if the 
fuselage is free to move, two additional equations are required to obtain 
the two additional unknown quantities cty and z r . These additional 
equations are those expressing the dynamic equilibrium of the forces 
on the fuselage, namely. 


+ Lt - m-£>^z r - Ax Oj.) = 0 
Mw - x t Lt - (r f 2 + (Ax^mfa,, + Ax xaf'z r - 0 




(80) 


where the wing lift and pitching moment can be expressed In terms of the 
lift and moment distributions l and m as 


or 



hr “ < lL@J| z | + 4®J{ W } 


L@J] z j + 





(81) 


Mw = <1 


(82) 
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where [_@)J > |_(§)J > j_(§)J ' and [_@)J 8X6 the rowS °" btained ^ 

postmultiplying the rows of equation ( 8 l) by the square and rectangular 
matrices of equation (66). In equation (80) the fuselage lift and moment 
have been neglected; they can easily be expressed in terms of and 

z r and included, if desired. 

The equations for the tail deformation Ax and those for the 
overall normal force and pitching moment can be combined with equa- 
tion (7 6) as follows: For the sake of def ini teness it will be assumed 

that n stations on the wing are considered, including the one at the 
root, so that there are 2 n + 1 unkn own quantities, and that in the 
column matrices defining applied loads, deformations, and so on, the 
values at the root of the wing are written at the top: 



The first and (n + l)th equations of the system defined by equa- 
tion (76) express only the trivial fact that the structural deformation 
at the wing root is zero. They are replaced by equations (82) and (79), 
which are adjoined to the system, to yield the combined equation of 
motion 



(83a) 
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If "the square matrix on the right side of this equation is designated 
by \k\> the rectangular matrix "by [b] , and the quasi-unit matrix on 
the left side of the equation by [l'J , the equation can also be written 
as 




(83b) 


The matrices 




are the same as the matrices 



to 


of equation ( 7 6 ), except that the first rows of the latter. 


which are all zero, are replaced as follows: 


First row of - 

Replacement 

Quantity added to the 
leading element 

© . 

First half of (l^J 

F x (k) + k 2 

id) 

© 

Second half of |_©J 

nu Ax „ 

f 2 00 - f 0 

m 

© 

First half of |^(S)J 

-xt Fide) - k 2 

ID 

© 

Second half of |^( 12 )J 

f 2 « + ^ 1. 2 

pfc) 

2 W 

© 

L©J 

0 

© 

L@J 

0 






MCA TN 3910 


69 


Also, the elements of the last rows and columns of the matrices £a] 
and [B] are zero, except for the following: 


A l,2n + 1 “ F 3 ( k > " ^ V k > 


^n+l,2n+l -x t 


K, 


F 5 (k) “ ^ *\( k )| 


* 211 + 1,1 " “ K 1 F l ( k ) 

*2n+l,n+l = “ K 1 F 2( k ) 


*2n+l, 2n+l K 1 


Ko 

F,(k) - F^k) 

Ki 


B l,rri-1 = F 5^ k ^ 
En+^n+l = -x t F 5( k ) 
^+1,^1 * “ K 1 F 5^ k ^ 


Solution of the Equations of Motion 


For the purpose of calculating the desired transfer functions, 
equation ( 83 b) may be solved directly for a given value of q as a set 
of linear algebraic equations with coefficients given by the matrix 


[l’] - q[A] and with "knowns" given by the column matrix q[B] 


|where [a], [b] , and 


)• 


are functions of kl. The result is a 


column matrix of the unkno wn amplitudes of the motions of the airplane. 
Xf this column is calculated for several values of k in the range of 
Interest, these amplitudes, considered as functions of k, are transfer 
functions from the gust to the motions. 


This column matrix can be substituted into equation (66) and the 


resulting column matrix 



substituted into equation ( 7 ^), to yield 


a column matrix of bending and twisting moments which again, considered 
as a function of k, represents transfer functions from the gust to these 
moments. A set of transfer functions for the vertical shear could be 
calculated similarly from the relation 



The stress at any point of the structure can be assumed to be given 
by a linear superposition of the bending moment, twisting moment, and 
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vertical shear at the given station,, if elementary beam theory is used. 
If elementary "beam theory cannot he used because of the interaction of 
bending and torsion stresses or because of shear lag, the stress at a 
given point can be expressed as a linear superposition of moments and 
shears at other stations as well as the given station. In either case, 
the transfer function for the given stress is then produced by the same 
linear superposition of the transfer functions for the corresponding 
moments and shears . 

It may be noted that, at zero frequency, solution of equation (83b) 
yields the static aeroelastic deformations and thus permits the calcu- 
lation of the changes in the lift distribution and the shift of the 
aerodynamic center that result from static aeroelastic aetion. Also, 
inasmuch as this equation completely describes the dynamic behavior of 
the airplane, the speeds at which aeroelastic instability phenomena 
occur can be calculated from it, although such calculations are beyond 
the scope of this paper. Suffice it to point out that for such a calcu- 
lation the degrees of freedom of the body must be eliminated first, as 
a result of the way in which the problem has been set up. This elimin- 
ation can be effected readily by considering the firBt and (n + l)th 
rows of [A], but with A 1L , A^ n+1 , Vi,! Vl,n+1 replaced 

by 0 . If these rows are premultiplied by 

A 11 A l,n+1. ~ 1 

An+1,1 A n+l,n+l 


and used as the first and (n + l)th rows of a matrix which is otherwise 
a unit matrix, and if this resulting matrix is referred to as £l"] > 
then equation (83b) can be written for this homogeneous case as 


[M [1"] - 

L ■ lAoj l J 


The products [l'] [l"] and [a] [l"] will now have two null rows and 
columns each, which correspond to z r and a^. If these rows and columns 


are deleted and z r and a„ are deleted in 4 


a L yielding a column 
Ax 


< a > , the remaining matrices are nonsingular, so that they can be inverted 


and the preceding equation can be written as 



ION 


NACA TN 3910 


71 


[w - ifrOO]] 


z 
a 
I Ax 


4 a ^ 4 0 > 


where [e] is the identity matrix, and 

[boo] 


m 


The horizontal braces designate the fact that the null rows and columns 
have been deleted. 

Equation (8i{-) is in the canonical form for the calculation of 
eigenvalues. If k is set equal to zero and the eigenvalues of jjD(0 )j 

are calculated by iteration, expansion of the determinant, or any other 
suitable method, the lowest real and positive one represents the value 
of the dynamic pressure at divergence. For swept wings the value lowest 
in absolute magnitude is usually negative and Is therefore of no practical 
significance, although it is often used as an index of the aeroelastic 
behavior of the airplane. 

This calculation can be repeated for various positive values of k, 
the first few eigenvalues being obtained for each. The results, which 
will generally be complex, can be plotted against k. When any of the 
eigenvalues becomes purely real, it represents a dynamic pressure at 
flutter, and the corresponding value of k represents the reduced fre- 
quency at flutter. (This statement is true only if the structural damping 
Is zero; such damping effects can easily be included, but the details of 
the process are beyond the scope of this paper.) 


DINAMICS OF THE LARGE FLEXIBLE AIRPLANE 


For the large flexible airplane the fundamental proposition of 
power spectral analysis, that the output power spectrum of a system is 
the product of the absolute square of the transfer function and the input 
power spectrum, is no longer valid if the Input is considered to be the 
gust intensity at a point. Nor can the output power spectrum be expressed 
directly in terms of an effective input spectrum, as in the case of the 
rigid airplane, where this simplification resulted from the assumption 
that the indicial-response function was expressible as the product of a 
function of time alone and a function of distance along the span alone. 

In the first section of this part of the paper the statistical 
problems involved in an analysis of the response of a large flexible 
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airplane are considered. The nature of the generalized transfer functions 
required for this purpose is described, and the means whereby they are 
combined with the input spectrum are indicated. 

The second section is concerned with an extension of the method 
outlined in the preceding part to the case of the large airplane. 

Although fundamentally the dynamic aspects of the problem are unchanged, 
and although the longitudinal and lateral degrees of freedom can still 
be separated, a direct application of the approaches outlined in the 
preceding parts of the paper to the large airplane requires consideration 
of the entire wing, rather than only one half of the wing. Little addi- 
tional computing time is then required to treat the lateral and longitu- 
dinal degrees of freedom simultaneously, and the necessity of combining 
the results of two separate analyses is obviated. However, attention 
can still be confined to one half of the wing by using the technique 
outlined in the discussion following equation (55) and in the section 
headed "Combination of the Results Obtained From Analyses of the Longitu- 
dinal and Lateral Degrees of Freedom. " Basically this technique consists 
in splitting the influence functions of concern into symmetric and anti- 
symmetric parts and using one part for an analysis involving the longitu- 
dinal degrees, and the other in an analysis involving the lateral degrees. 
If this approach is adopted, separate analysis of longitudinal and lateral 
degrees of freedom is still preferable. 


Extension of the Statistical Approach 

The power spectrum and, hence, the mean-square values of the 
responses such as the stress at a given point on the wing of a large 
flexible airplane due to flight through turbulent air can be calculated 
in several ways. Perhaps the most direct of these consists in using the 
basic approach outlined in the first part of this paper and starting 
with an expression for the instantaneous value of the stress in temns 
of a suitable Indicial-response influence function, namely, 

p 00 pb/2 

a(t) =J J ^('tijy) v(u(t-t L },y) dy dt x (85) 

The function h^(t,y) Is, as before, an indicial-response Influence 

function or Green's function for the partial differential equation. It 
relates the stress as a function of space and time coordinates to the 
applied loads,, which are also functions of space and time coordinates. 

The essential difference between the' problem considered in this section 
and those treated in the first two parts of the paper is that this func- 
tion can no longer be expressed as a product of a function of t alone 
and y alone. In some cases - In a modal approach, for instance - it 



NACA TN 3910 


73 


may be expressible as a sum of several such functions, and then the 
approach used previously for the bending of a large rigid airplane free 
to move vertically may be adopted, as has been done in reference 13. 
However, in this section the case is considered in which even this 
simplification cannot be made. 

The correlation function for cr can be calculated directly from 
equation (85), and for the case of axisymmetric turbulence it is 


V t) 2 ^(V y l) *£(V y 2) ** (j u2 ( T+t i- t 2) 2+ (y 2 - y i} 2 ) ^2 dt i dt 2 W) 

Hence, the power spectrum of a can be obtained by calculating the 
Fourier transform of this function and is 

pb/2 rib/2 ^ 

9ff(a5) J b / 2 J b / 2 ^ (®' y l) ^(^ y 2) ^w(^ |y2- y lf) ^1 ^2 ^ 

where ^(03,11) is the two-dimensional spectrum considered previously, 

H^(ca,y) is the Fourier transform of h^(t,y) with respect to time, 
and the asterisk designates a complex conjugate, as before. Thus, the 

function H^(cn,y) is a transfer function from sinusoidal vertical gusts 

(of width dy) impinging on the wing at a given station y on the wing 
to the stress a, or a Green's function for the ordinary differential 
equation (with the quantity cu as a parameter) that relates the stress 
amplitude as a function of the space coordinate y to the amplitude of 
the applied sinusoidal gusts . 

The term H a (a>,y-|J H^o^y 2 ) 1 - n equation (87) is complexj however, 

the imaginary part can be ignored because it contributes nothing to the 
integral as a result of the fact that cfJ^ depends only on |y 2 

In a manner analogous to the one employed in the first part of this 
paper, the double integral in equation (87) can be evaluated by intro- 
ducing the function H^(o>, T}) which takes the place of the function r(^) 

used previously and is defined by an ■ autoconvolution of H^(o3,q): 
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7^ 

■where designates the real part. Hence, 

cp (oi) = f $*(<»> n) dT i ( 8Q ) 

^ o 

Another approach consists in using the double Fourier transform. 
$ W (A) of the axisymmetric input correlation defined in the first part 

of this paper. In terms of this function the correlation function t ff ( T ) 

can be written as 


t;( T > " Ijz lj / z /_ h o( t:L ’ y ^ e 1 ^ lU ( T+t i- t 2) +A 2 ( y 2‘ y 0] t v (js 1 z +\ s s j ax x dx 2 ay x dyg dtj. atg 

‘ iflfl ° lhUT 


where the function 


^(^1^2) = / ro / b / 2 h 0^> y) e" 1 ^^) ay dt 


represents the Fourier transform (with respect to y) of the function 
H^(cD,y), that is, 

' 1/2 _i 


= f °^ 2 e -1 * 27 ^(H u >y) ay 


(89) 


Hence, 



With the approach outlined in the next section, the function 
H cy (cn,y} can be calculated either directly or indirectly, by first calcu- 
lating the transfer function from that gush to the lift distribution and 
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then the transfer function from the lift distribution to the stress. 
For the indirect method , 



where the function Hj/cd) G(rj,y) is the influence function for the local 

lift or the Green's function for the aerodynamic problem involving s in us - 
oidal gusts considered in the first part of this paper; the symbols T] 
and y in G(T],y) are interchanged, however, so that the function now 
defines the contribution of a gust at station y to the lift at sta- 
tion T]. The transfer function relates the (sinusoidal) stress 

at the given point to a unit concentrated (sinusoidal) normal force 
acting at station tj. 

With this indirect method, the power spectrum for the stress at a 
given point can be calculated by starting with the power Bpectrum for 
the lift distribution calculated in the first part of the paper. For 
this approach o(t) may be written as 

p°° pb/2 7 

° (t) MV y ) ‘(‘-V) ay at i 


so that 


*a(r) 


r>°° pco To/2 p-b/Z 2 

-«J-»J.b/2J-b/2 M T+t l-W 2 ) ^2 dt l dt 2 


and 


p b/2 rb/2 . 

^ ff (“) J b y 2 J b y 2 ^ (“^l) <^1 ^2 ( 9 !) 


>b/2 

-b/2 


where 


The Fourier transform of this correlation function is, then. 



( 92 ) 
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(See the section headed "Generalized Aerodynamic Influence Functions in 
Unsteady Flow. 11 ) Hence, 




,2 pl>/ 2 pb/2 * 

H * “)' J t / 2 J Tj/g ( m ' y i) M®' y a) ^2 


(95) 


If cd (cc) is given by equation (35)> the double integral can be 
• "€ “ 
expressed in terms of single integrals as follows: 

qp a (o>) = |h^(o>)| 2 [(1 - K) 2 |G(a3)| 2 cp We (a>) + 2K(l - K)R |g(cd) G'*(a>)J + 


K c 


r b 

■ G"(cu, rj) cp v (<m,Ti) dTj 


(94) 


where 


,b/2 


, H-(co,y) y(y) dy 
-b/2 

1 

G'(ai) = J 


b/2 

H£(<n,y) cp w _*(a>,y) 7 ( 7 ) <3y 

■b/2 e 


as before, R {} 
imaginary part. 


»"(«) s 2 |h |Hj(co,y)J R|Hp(o),y+ri)| + 

l|Hj(<n,y)| l|Hj(o),y+Ti)|J /(y) r(y+Ti) dy 

:{} 


designates the real part, and 1^ f designates the 


Equations (88), (90), and (93) thus represent three methods of 
obtaining the power spectrum of the given stress. One requires a trans- 
fer function from the local gust intensities directly to the stress, an 
autoconvolution of this transfer function, and the spectrum of turbulence 
define! by equation (l6)$ the second requires a two-dimensional spectrum 
of the turbulence defined by equations (l7 a ) and (17b) and a Fourier 
transform of the aforementioned transfer function with respect to yj 
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the t hi rd utilizes an autoconvolution of the transfer functions from 
local concentrated loads to the stress and the spectrum, for local lifts 
due to vertical gusts calculated in the first part of the paper. The 
choice of approach depends to some extent on the information available, 
but is largely a matter of individual preference. 

The analysis in this section has been based on the premise that 
both halves of the wing would be treated simultaneously. Attention can 
be confined to one wing by using only the symmetric parts of the influ- 
ence functions in an analysis involving the longitudinal degrees of freedom, 
anrl only the antisymmetric parts in an analysis involving the lateral 
degrees of freedom. The symmetric part of an influence function for a 
unit concentrated load or gust acting at station y is the response func- 
tion for two loads or gusts of l /2 unit intensity acting at stations y 
and -y, respectively. Similarly, the antisymmetric part is the response 
for a load or gust of l /2 unit intensity acting at station y and an equal 
and opposite load or gust acting at station -y. If this approach is 
chosen, seme of the integral expressions given in this section assume 
slightly simpler forms. For instance, equation ( 89 ) can then be written 
for the symmetric part as 



cos Agy ^(X x U,y) dy 


and for the antisymmetric part as 



sin Agy H^A L U,y) dy 


Calculation of the Required Transfer Functions 

Depe nd ing on which of the methods outlined in the preceding section 
is used, one of two types of transfer functions is required - either the 
one from local gusts to the stress of interest, or the one from the local 
lift to that stress. Both of these functions differ in several respects 
from those considered in the preceding part of this paper. 

For either type of transfer function the structural deformations of 
the wing under concentrated loads as well as under distributed loads 
■will be required, so that the numerical-integration schemes used in the 
preceding part have to be modified to seme extent. Also, it may now be 
preferable to consider the entire wing (rather than the one semispan 
considered in the case of the small flexible airplane) , so that the 
various transfer functions are asymmetric. With the degrees of freedom 
considered in the analysis thus doubled, it becomes preferable to include 
one additional degree of freedom, roll, rather than to perform two sepa- 
rate analyses for symmetric and antisymmetric motions with, respectively. 
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one and two less degrees of freedom, and then to combine the results. How- 
ever, if the alternative approach of splitting up the influence functions 
into symmetric and antisymmetric parts for use in two separate analyses is 
adopted, only one semispan need be considered, and the results presented 
in this section can then be simplified to a large extent. 

Before discussing the modifications required to extend the dynamic 
analysis outlined in the preceding part of the paper to the large air- 
plane, it might be pointed out that chordwise deformations (deformations 
parallel to the chord) will again be ignored, ._ Again, they can readily 
be included by a straightforward extension of the approach used here if 
it is felt that they may be significant in any given case. If they are 
included, however, yawing and possibly also sideslipping motions can 
probably no longer be ignored, because they may give rise to large forces 
in the chordwise direction. If the entire wing is treated, these two 
additional degrees of freedom can readily be included, and all longitu- 
dinal and lateral degrees of freedom are then treated simultaneous lyj 
if two separate analyses are performed for the longitudinal and lateral 
degrees of freedom, the symmetric and antisymmetric parts of the influ- 
ence functions being used and only one sepiispan being treated, these 
additional degrees of freedom enter only into the lateral analysis. 

The structural deformation due to local (concentrated) loads can 
be obtained in several ways. If measured influence coefficients are 
used, they pertain precisely to such loads and, in fact, must be modified 
before they can be used for distributed loads (see ref. 12) so that it 
is necessary only to use the unmodified coefficients. 

If the deformations are to be calculated in a manner similar to that 
employed in the preceding part, the integrating matrices must be replaced 
as follows: 


|W 



where 


H 

II 

(q. > p) 

HlOJ 

II 

o* 

(q. = p) 

•cT* 

!^ 

II 

o 

(q. < p) 


and 


(I) 2 [ll] by |[ll'] 
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where 


Upq' - yq - yp 


(q. ^ p) 


Ilpq' =0 (q < p) 

The factor of 1/2 for Ipp ' constitutes an approximation which implies 
fairing through a discontinuity . If this approximation is to be avoided, 
the deflections due to unit concentrated loads (the structural influence 
coefficients) can be calculated directly from simple beam theory, in 
which case the limits of integration take care of the discontinuities. 
Thus, for instance, for an unsvept wing, the normal deflection and twist 
at yp due to a unit concentrated load and torque, respectively, at y^ 
are 


Z P<1 

z pq. 

c'pq 

°pq 


„ r y Pf yt yq - y 


r p r 

Jo Jo 


dy dy' 


El(y) 

^ f y ’ yq 7 y ay ^ 


= r q r ■ 

Jo Jo 

_ r p 1 

J 0 C3J(y) 


Ei(y) 


dy 


-I, 


y q 


QJ(y) 


dy 


+ 



yg - y 

El(y) 


dy 


(y p £ yq} 

(y P > yq) 

(y P £ yq) 

(yq < yp) 


The concentrated loads under consideration arise as follows: For 

the transfer functions relating local lifts to the desired stress, the 
local lifts may be considered to be concentrated loads of unit magnitude, 
associated with concentrated torques of magnitude epC. Equation (83b) 
can then be written as (see also eq. ( 76 )) 




(95) 


where Is a diagonal matrix of the values of w defined by equa- 
tion ( 63 ), and where the matrix represents either the four 

influence-coefficient matrices for z and a due to concentrated loads 
and torques, or the square matrix of equation ( 75 ) with modified inte- 
grating matrices, as discussed in the preceding paragraphs. 
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It may be noted that equation (95) uow represents not onfe set of 
simultaneous equations hut several, all having the same coefficients 
but different sets of knowns (as defined by the columns of the matrix 

on the right side) and, hence, different sets of unknowns (the columns 


of the matrix 



This situation is. due to the_ fact that the func- 


tions under consideration are, in effect, the responses of the airplane 
as a whole to sinusoidally varying concentrated loads and are different 
for each location of the applied load. 


Once this equation has been modified to take into account the 
overall body motions and tail deflections (as explained in the preceding 
part) as well as the rolling motions (as explained in the following 


paragraphs), it can be solved to yield the unknown values of 


z 

a . 


From 


LA*] 

these values the bending and twisting moments, as well as the vertical 
shears, can be calculated and added to those due to the concentrated 
loads. When combined linearly as required for the desired stress, these 


moments and shears yield the desired transfer functions 


Hj(cu,y) . 


If the transfer function directly from the local gusts to the desired 
stress is to be determined, the response of the airplane to the lift 
distribution induced by a sinusoidal gust of width dy acting at sta- 
tion y must be calculated. This lift distribution is the Green's func- 
tion considered previously. If it is represented by the relation 

E^(cu) G(y,q), with G(y, t}) defined by the approximation given in equa- 
tion (34), the concentrated loads arise from the delta function in that 
expression. The right side of equation (95) becomes, In that case, 


4(1 - K) H^(co) [>] 


|V(y) 7 (ti)1 


1( e i c |jV(y) rOi)] 


R + qbK H^(cd) [r'] 



where jffj represents the square matrix of equation (75), and [r'3 

the one discussed in connection with equation (95) • Again, several sets 
of simultaneous equations are implied. Their solution (after modifica- 
tion for overall motions and tail deflections) yields values of z and 

a from which the transfer functions H^(o>,y) can be calculated. 
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The extension of the method of the preceding part to the calculation 
of the deformations on both wings is straightforward. Essentially, dis- 
tributed lifts and torques now have to be calculated for both wings and 
integrated with matrices which can be assembled from those used for one 
wing alone. No new problems arise in this process, so that it need not 
be discussed further. 

The. inclusion of rolling motion, however, is not so straightforward. 
One method consists in replacing all values of z in equation (83a) or 
its equivalent by z + 0y and then reducing the columns involving this 
quantity by the following relation (which assumes that the new unknown 
quantity, the roll angle 0, is listed at the end of the column): 



1 

I 

1 

I 

' 


1 


°; < 
1 
1 

y 



z 

1 

! 

I 

1 

> 









0 



a 

l 




Ax 


1 



0 





L 


An additional equation must then be joined to the set, 
equation of equilibrium in roll 


namely, the 



Ijr0 - Mp^.0 = 0 


(96) 


where I x is the inertia in roll of the fuselage and empe nnag e alone 
inasmuch as the inertia effects of the wing are included in Z(y), and 
is the coefficient of damping in roll for the empennage. For most 
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cases both of these contributions are negligible. If they are ignored, 
equation (96) can be written in matrix notation as 

Lnj(i) - 0 

where [nj is now a matrix which serves to perform the integration 

required in equation (96). This condition can then be adjoined to the 
other equations of the set in the same manner as equations (82) were 
adjoined to the set in the preceding part. 

The result, again, is a set of simultaneous equations for z, a, 

Aa, and 0 from the solution of which the desired transfer functions 
can be obtained as outlined in the preceding paragraphs. Also, as before, 
once the unknowns z r , a^, and 0 are eliminated from the set, the 
divergence and flutter speeds can be calculated by conventional matrix 
operations; these speeds will then pertain to an airplane free to move 
vertically as well as in pitch and roll and, hence, will include divergence 
and flutter speeds in antisymmetric as well as symmetric modes. 


DISCUSSION 


Some Implications of the Assumptions Concerning the Nature 
of Atmospheric Turbulence 

Hie turbulence was assumed to be homogeneous in order to make the 
problem stationary in the statistical sense and thus permit the use of 
the mathematical techniques developed for such problems. In a practical 
sense, turbulence can be homogeneous only in a limited body of air. The 
assumption thus implies that the dimension of this body of air along the 
flight path is large compared with the distance traversed in the reaction 
time of the airplane, which in the case of load studies is of the order 
of the time to damp to one-half amplitude, but in the case of motion 
studies may be much larger. Obviously, the greater the body of air, the 
greater the reliability with which the loads and motions can be predicted 
(in a statistical sense) for one run through it. In general, turbulence 
at very low altitudes, which may be influenced significantly by the con- 
figuration of the ground, and turbulence in thunderstorms, may not be 
sufficiently homogeneous for the purpose of this type of analysis, but 
other types of turbulence are likely to be substantially homogeneous over 
sufficiently large distances. 

Isotropy was assumed in order to permit the required two-dimensional 
correlation functions to be expressed simply in terms of the one-dimensional 
correlation functions. For sufficiently short wave lengths all turbulence 
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is isotropic, but for long wave lengths it can be isotropic only if it is 
homogeneous (both in the plane of the flight path and perpendicular to 
it) . The condition of axisymmetry, which is sufficient for most of the 
results presented herein, is less restrictive than isotropy inasmuch as 
it does not specify the variation of the characteristics of the turbulence 
in the vertical direction. In practical problems, if the turbulence may 
be assumed to be homogeneous, the Conditions of axisymmetry and isotropy 
are likely to be satisfied to a sufficient extent to permit the use of the 
approach presented herein for all but very long wave lengths. The wave 
length at which it ceases to be valid depends on the size of the body of 
air under consideration, being larger for a large body. 

Taylor's hypothesis (to the effect that a space displacement Ax 
along the flight path may be identified with a time displacement t = Ax/U 
in the gust correlation functions) implies that the variation in gust 
intensity that prevails along the flight path at any instant will remain 
substantially the same until the airplane has traversed the given body of 
air. The required correlation functions for atmospheric turbulence are 
thus in the nature of space correlation functions (rather than time corre- 
lation functions) and have been considered as such. The statistical char- 
acteristics of the turbulence are then independent of the speed at which 
it is traversed. Clearly, whether or not this hypothesis is valid depends 
on the flying speed of the airplane. On the basis of present knowledge 
no definite lower limiting speed can be quoted. However, indications are 
that the hypothesis is valid for flying speeds greater than about 100 or 
200 feet per second. The effect of finite flying speed on the gust corre- 
lation function can be expected to be most pronounced for large distances, 
where the correlation is weak, so that the effect on the various spectra 
is likely to be small and to occur at the longest wave lengths, where, as 
previously mentioned, the spectrum is somewhat uncertain for other reasons 
as well; this effect is thus more likely to be significant for large than 
for small airplanes. 

The particular correlation function used herein for the calculations 
of the "averaged" correlation functions and spectra has certain theoretical 
shortcomings - primarily that the associated spectrum does not decrease 
rapidly enough for very short wave lengths. However, it does appear to be 
adequate to represent the available information concerning the spectra of 
atmospheric turbulence (see ref . 17 , for instance), because the behavior 
at very short wave lengths is relatively unimportant , inasmuch as airplanes 
cannot respond to them, and the behavior at very long wave lengths is 
usually in doubt by reason of the nonhomogeneity of actual turbulence. 3h 
the intermediate range of wave lengths, this correlation function appears 
to be quite satisfactory. 

The parameter L* (the integral scale of turbulence) used herein is, 
for practical purposes, a largely fictitious quantity, inasmuch as it is 
proportional to the values of the gust spectrum for infinite wave lengths. 
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which, in view of the uncertainties in the values of the spectra at large 
wave lengths, have little physical significance. Therefore, at present, 
insufficient information is available to give a value for L* to be used 
in connection with the numerical results calculated herein, although a 
value of 1,000 to 2,000 feet appears to be appropriate. As more informa- 
tion concerning the spectrum of atmospheric turbulence becomes available, 
more definite values can be deduced by fitting an analytical expression 
of the type used here to measured results in the range of frequencies of 
primary interest, and then using this expression as a means of obtaining 
a value of L* by extrapolation of the measured results to infinite wave 
lengths (zero frequency). 


Considerations Pertinent to the Application of Stationary- 
Random-Process Tec hni ques to Gust-Load Problems 

The purpose of this section is to point out how, in principle, the 
results of analyses of the type outlined herein may be used in overall 
load analysis and, hence, in the design of an airplane. 

Consideration is confined in this paper primarily to the power spectra 
of the motions and stresses of interest. As pointed out in references 1, 

2, 3, and 11, for instance, a great deal of statistical information of 
direct interest can be obtained from the power spectrum. For instance, if 
the random process of concern (say, the given stress as a function of time) 
has a Gaussian probability distribution, the expected number of peaks at 
or beyond a given level in a given period, of time can be calculated very 
simply from the integral of the spectrum and its second and fourth moments. 

The results obtained in this manner pertain to continued flight in a 
given body of turbulent air. They have to be generalized by determining 
the likelihood of flying through turbulence of the given characteristics. 
(See ref. 18.) The probability of exceeding a given stress level during 
the expected life of the airplane while flying through atmospheric turbu- 
lence can thus be calculated in straightforward fashion. To this proba- 
bility must then be added the probability of exceeding this level in 
maneuvers, landings, and, possibly, also in turbulence due to thunder- 
storms, because in view of the possibly nonhomogeneous character of turbu- 
lence in thunderstorms and the possibly nonlinear nature of the aerodynamic 
forces incurred while flying through them, the techniques" used herein may 
not be applicable to flight through thunderstorms, and a separate analysis 
may have to be performed. 

Although the available information concerning atmospheric turbulence 
is inadequate to permit of any definite conclusion, the results obtainable 
with the approach outlined herein may turn out to be most significant for 
the prediction of the low- and medium-amplitude stress cycles which are 
important to fatigue studies; their validity for or contribution to the 
prediction of very severe loads remains to be seen. 
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CONCLUDING REMARKS 


The statistical approach to the problem of calculating the dynamic 
responses and the stresses of an airplane subjected to continuous random 
atmospheric turbulence has been extended in several respects; basically, 
only the assumptions of linearity, that is, of small motions and deforma- 
tions, as well as homogeneity and axisymmetry of the turbulence are 
retained. 

The first problem considered was the effect of spanwise variations 
of the instantaneous turbulent velocities on the lift and moments due to 
turbulence. The mean-square lift has been shown, to be reduced consider- 
ably if the span of the airplane is relatively large compared with the 
integral scale of turbulence. The shape of the spectrum of this lift is 
affected relatively little by spanwise variations of gust intensity, 
except at very high frequencies, if the decrease in the effective mean- 
square intensity is taken into account. The effect of sweep on the mean- 
square lift and its spectrum has been shown to be small for wings with 
a given distance fhom root to tip. 

If the variation of the instantaneous velocities is taken into 
account, the rolling moment to which the airplane is subjected can be 
calculated. The mean-square rolling moment has been shown to be propor- 
tional to the ratio of the wing span to the integral scale of turbulence 
for small values of that ratio. Similarly, expressions for the mean- 
square values and the power spectra of the local lift, the bending moments, 
and the pitching moment have been given. For some of these forces the 
required aerodynamic information cannot be calculated by existing methods. 
Therefore, certain approximations, based on experience with steady aero- 
dynamic forces and available knowledge concerning unsteady forces, had to 
be made for the aerodynamic influence functions in unsteady flow. 

The next problem considered was the dynamic response of a rigid air- 
plane to random turbulence . This problem had previously been treated 
for the case of an airplane free to move only in the vertical direction 
and small enough so that variation of the turbulent velocities along the 
span could be neglected. In the present paper the response of an airplane 
in three longitudinal degrees of freedom was considered; calculations were 
made which suggest that the inclusion of deviations from the mean hori- 
zontal motion is superfluous in gust-load calculations. For the remaining 
two longitudinal degrees of freedom, the mean-square normal and angular 
acceleration have been shown to be functions of only two parameters other 
than the mass ratio and scale parameter of the single-degree-of -freedom 
case, namely, dimensionless forms of the short-period frequency and of 
the time to damp to one-half amplitude. An indication is given of the 
manner in which the results obtained in connection with the first problem 
can be used to extend this dynamic analysis to the case in which variations 
of the turbulent velocity along the span have to be taken into account. 
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The last problem treated was the dynamic response of a flexible 
airplane, including vertical motion, pitch, and, when necessary (as 
when spanwise variations in gust intensity are taken into account), roll. 
Horizontal and lateral (yawing and sideslipping) motions were disregarded 
because they do not generally affect the wing stresses due to vertical 
gusts. A method which represents an extension to the dynamic case of a 
numerical-integration approach to the static aeroelastic problem has been 
outlined for the analysis of the problem at hand. The modifications 
required in the basic statistical approach and in this method of dynamic 
analysis in order to treat the case in which spanwise variations of the 
gust intensity are important have been discussed. 

Although most of this analysis has been confined to the vertical 
component of turbulence, it has been Bhown that the simultaneous action 
of longitudinal, vertical, and lateral gusts on the wing stresses (with 
due allowance for the fact that vertical gusts affect both the longitudi- 
nal and the lateral motions of the airplane) can be taken into account by 
simply adding the power spectra of the various contributions, provided 
the turbulence is isotropic; the cross correlations or spectra have been 
shown to vanish either by the symmetry or antisymmetry of the influence 
functions involved or as a result of the statistical independence of 
mutually perpendicular velocity components. 

The approach presented herein thus furnishes a foundation for the 
prediction of the statistical properties of the stress experience of a 
given airplane once the appropriate statistical characteristics of the 
atmosphere have been determined. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 5, 1956. 
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TABLE 1.- DEFINITIONS AND NUMERICAL VALUES OF 
STABILITY DERIVATIVES 


Derivative 

Definition 

Numerical value 
for the example 


-*§(% + Cd ) 

-l.lt-30 

ZxL 

+ C l) 

-0.0955 

Xw 

f ( C L - °»a) 

0.0016 

X u 

■^(°Du + °d) 

-0.0097 


qSc p 

-0.0235 


qSc 2 r 
2IyyU 2 ^ 

-0.0013 

M u 

PK * ^ 

0 

Mq 

q.sc 2 

2IyyU % 

-1.920 


W, lb 30,500 

U, fps 660 

c, ft 10 

Altitude, ft . 20,000 
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Figure 3 • - The power spectrum of the averaged vertical component of 

turbulence . 


















